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Under appropriate conditions, mixtures of cationic and neutral lipids and DNA in water condense 
into complexes in which DNA strands form local 2D smectic lattices intercalated between lipid bi- 
layer membranes in a lamellar stack. These lamellar DNA-cationic-lipid complexes can in principle 
exhibit a variety of equilibrium phases, including a columnar phase in which parallel DNA strands 
from a 2D lattice, a nematic lamellar phase in which DNA strands align along a common direc- 
tion but exhibit no long-range positional order, and a possible new intermediate phase, the sliding 
columnar (SC) phase, characterized by a vanishing shear modulus for relative displacement of DNA 
lattices but a nonvanishing modulus for compressing these lattices. We develop a model capable of 
describing all phases and transitions among them and use it to calculate structural properties of 
the sliding columnar phase. We calculate displacement and density correlation functions and x-ray 
scattering intensities in this phase and show, in particular, that density correlations within a layer 
have an unusual exp(— const. In^ r) dependence on separation r. We investigate the stability of the 
SC phase with respect to shear couplings leading to the columnar phase and dislocation unbinding 
leading to the lamellar nematic phase. For models with interactions only between nearest neighbor 
planes, we conclude that the SC phase is not thermodynamically stable. Correlation functions in 
the nematic lamellar phase, however, exhibit SC behavior over a range of length scales. 
Pacs; 61.30.Cz,61.30.Jf,64. 70.Md 
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I. INTRODUCTION 

The search for non- viral vectors for transport of DNA 
across cell and nuclear membranes in gene therapy has 
led to the study of DNA-cationic-lipid complexes [Q. 
DNA and mixtures of neutral (zwitterionic) and cationic 
lipids dispersed in water self-assemble into spheroidal 
complexes that can attain micron sizes near the isoelec- 
tric point where there is compensation between DNA 
and lipid charge. Fluorescence tagging of DNA and 
lipids indicate that both species are dispersed uni- 
formly throughout the complexes. X-ray scattering ex- 
periments reveal two bulk structures for the com- 
plexes with close association between DNA and lipids. 
When helper lipids that favor spontaneous curvature of 
lipid membranes or that decrease membrane charge den- 
sity are added, a hexagonal inverted micellar structure 
forms with DNA molecules captured in the water holes of 
the hexagonal lattice Q . When only neutral and cationic 
lipids are used, DNA is intercalated in galleries between 
lamellae of a lamellar lyotropic phase formed by the lipids 
l^-Q as depicted schematically in Fig. |l[ This article 
will investigate the possible equilibrium phases of these 
lamellar DNA lipid complexes. It will focus primarily on 
the properties of one phase, the sliding columnar phase 
1^,0 , characterized by strong orientational but weak po- 
sitional correlation between DNA strands in neighboring 
galleries. 




FIG. 1. Schematic representation 

of lamellar DNA-cationic-lipid complexes. Parallel strands 
of DNA form 2D smectic lattices with lattice spacing d in 
galleries between flat lipid bilayer membranes with spacing a. 
DNA strands are aligned parallel to the x axis and the y axis 
is normal to the lipid planes. The height of the midpoint of 
the nth membrane above a point r = {x, z) in the a;2:-plane 
is -H"Scm(r). The nth DNA gallery with height -H'5na(i') lies 
between the nth and the n -|- 1st membrane [See. Eqs. (2.1) 
and (O)]. 



DNA molecules are semi-flexible polymers that when 
confined to a two-dimensional plane tend to form locally 
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aligned structures with a preferred inter-molecular sepa- 
ration that can be modeled as two-dimensional smectic 
liquid crystals J^-jlOt. Thus, lamellar DNA-lipid com- 
plexes can be viewed as a 3D stack of 2D smectics as 
depicted in Fig. |l|. This figure establishes the coordi- 
nate conventions that will be used throughout this paper. 
DNA strands align on average along the x direction, and 
the normal to lamellae is along the y axis. The normal 
to the 2D smectic lattices is along the z direction. The 
average spacing between lamellae is a, and that between 
DNA strands is d. The wavenumbers associated with 
these lengths are, respectively, fco — 2'k ja and go = 27r/d. 

If the lamellae are assumed not to be disrupted by dis- 
locations or rips, then the following possible equilibrium 
phases are easily identified: 

• Columnar (C) Phase: In this phase, DNA 
strands are aligned on average along the x axis, 
and their centers occupy positions on a 2D crys- 
tal lattice in the yz plane. Since the standard 
Coulomb repulsion between DNA strands favors 
staggering of smectic lattices in neighboring gal- 
leries, the columnar lattice is normally expected 
to be centered rectangular as observed in exper- 
iments on complexes in which membranes are in 
the L^' phase rather than the more disordered 
phase Q. We will, however, consider both simple 
rectangular and centered rectangular columnar lat- 
tices. A simple rectangular lattice may occur if the 
effective interaction between DNA strands in di/- 
ferent galleries is attractive, as is the case between 
DNA strands in solutions with polyvalent salts . 
The Columnar phase is favored at low tempera- 
ture. It has the same symmetry as a columnar 
discotic liquid crystal phase ||l2|. It is character- 
ized by a 2D elastic energy with a nonvanishing 
shear modulus for relative displacements of DNA 
lattices in different galleries and nonvanishing mod- 
uli for compression of both lamellae and DNA smec- 
tic lattices. It has long-range positional order in the 
2D yz-plane and associated Bragg peaks in its x- 
ray scattering profile at reciprocal lattices vectors 
Grm,m2 = [0, mifco/2, m2(7o] = {^.1,1^2) where mi 
and 1712 and integers and mi is even for a simple 
rectangular lattice and mi + m2 is even for a cen- 
tered rectangular lattice as shown Fig. ||. 

• Nematic Lamellar (NL) Phase: In this phase, 
the periodic positional order of the columnar phase 
is destroyed by dislocations in the DNA smectic 
lattices, but the long-range orientational order of 
DNA strands is maintained. This phase is char- 
acterized by a long- wavelength elastic energy with 
a lamellar compression modulus, an anisotropic 
lamellar bending modulus, and orientational rigidi- 
ties (Frank elastic constants) opposing spatially 
dependent variation of DNA alignment direction. 
Both the shear modulus for relative displacement 



of DNA lattices in different galleries and the com- 
pression modulus for DNA lattices vanish, and 
there is exponential decay of DNA positional cor- 
relations. The x-ray scattering profile of the NL 
phase exhibits lamellar power-law (2Z,0) peaks at 
Gl2i.o = (0, ?fco,0). If columnar-phase positional 
correlations are well developed, it will also exhibit 
Lorentzian peaks at Gmi,™.^ for m2 ^ as depicted 
in Fig. ||. If thermal fluctuations are sufficiently 
strong that these correlations are not well devel- 
oped in centered rectangular systems, then the x- 
ray scattering profile could exhibit a broad (0, 1) 
peak in the vicinity of (0, 0, go) rather than the ex- 
pected pair of (1, 1) and (—1, 1) peaks. In this case, 
the positions of the x-ray scattering peaks would 
appear to indicate a tendency to form a simple 
rectangular rather than the ground-state centered 
rectangular structure. 

• Isotropic Lamellar (IL) Phase: In this phase, 
orientational as well as positional order of DNA lat- 
tices is lost. Macroscopically this phase is identical 
to a multi-component isotropic lamellar phase. It 
is included for completeness and will not be consid- 
ered further in this paper. 

• Decoupled 2D Smectic (DS) Phase: This 
phase occurs only if there are absolutely no interac- 
tions between DNA lattices in neighboring galleries. 
Its elasticity and correlations are thus those of in- 
dependent 2D smectic lattices. Since there are al- 
ways interactions between galleries, this phase will 
not exist in real systems. It is, however, a useful 
limit to consider since systems with weak coupling 
between layers will behave as though they are de- 
coupled at sufficiently short length scales. 

In addition to the above phases, which are straightfor- 
ward to identify, there is the possibility of another phase 
with unusual properties: 

• Sliding Columnar (SC) Phase §,§: This phase 
has properties intermediate between those of the 
columnar and lamellar nematic phases. Its elastic 
energy is distinguished from that of the NL phase 
by the presence of a nonvanishing modulus for com- 
pression of DNA lattices. In-plane smectic correla- 
tions die off as exp(~const. In^ r) as a function of 
separation r rather than exponentially as in the 
NL phase. Correlations between smectic lattices in 
different galleries die off exponentially with layer- 
number difference. The x-ray structure factor ex- 
hibits power-law lamellar (0, Zfco,0) peaks and well 
defined DNA (mi, m2) peaks with m2 7^ that are 
sharper than the corresponding Lorentzian peaks 
in the NL phase, as depicted in Fig. ^. As in the 
NL phase, (1,1) and (—1,1) peaks may merge to 
produce a single (0, 1) DNA peak if thermal fluctu- 
ations are sufficiently strong. 



2 






(6,0) 
(4,0) 
(2,0) 

(-2,0) 
(-4,0) 
(-6,0) 



n qy 

O 

o 



27C/J 



•(5,1) 
■^•(3,1) 
•(1,1) 



^ Ill/a 



•(-1,1) 

•(-3,1) 
•(-5,1) 



n qy 






1 
• 




—U 




1 




• 

1 



(a) (b) (c) 

FIG. 2. Schematic representation of structure (top) and x-ray scattering profiles (bottom) of (a) the Columnar phase, (b) the 
Sliding Columnar phase, and (c) the Nematic Lamellar Phase. These figures assume that the preferred low-temperature phase 
is the centered-rectangular columnar phase. Modifications for a rectangular columnar phase are obvious. The indexing scheme 
for the x-ray intensities is the standard one for centered rectangular lattice. There are Bragg peaks in the columnar phase 
at G = (0, mifco/2, m2q'o) [denoted in the figure as (mi, 7712), where mi -I- m2 is even]. In the columnar phase, DNA strands 
form a 2D lattice with true long-range order and associated x-ray Bragg peaks. In the sliding columnar and nematic lamellar 
phases, there are power-law quasi-Bragg peaks at (2mi,0) and diffuse peaks at (mi,m2) for mi 7^ 0. Dislocations destroy 2D 
smectic-like order in the nematic lamellar phase to produce more diffuse (m\,m2) peaks than in the sliding columnar phase. If 
interlayer couplings are weak and thermal fluctuations sufficiently strong, the (±m, 1) peaks may merge to produce a diffuse 
peak centered at the rectangular-lattice (0, 1) position in the sliding columnar and nematic lamellar phases. The labelling of 
peaks in this figure, which follows conventional (x,y,z) ordering for the our choice of axes, is the inverse of that of Ref. Q. 
Thus (mi,m2) [e.g., (±1,1)] here corresponds to (m2,mi) [e.g., (1,±1)] in Ref. Q. 



Rather stringent conditions must be met before the SC 
phase can be thermodynamically stable. Thermal fluc- 
tuations must be strong enough to destroy the interlayer 
shear coupling present in the columnar phase but not so 
strong that dislocation proliferation destroys the smec- 
tic compressibility to create a nematic lamellar phase. 
The SC phase is stable only for temperatures T lying 
above a decoupling temperature Td at which the C phase 
becomes unstable and below a KT-melting temperature 
Tkt above which the NL phase becomes stable. Thus 
a necessary condition for the SC phase to be stable is 
Td < Tkt- In Sec. we will show that this condition is 
violated for the nearest neighbor models we consider here. 
Elsewhere we show that appropriately chosen in- 

teractions between further neighbor planes can stabilize 
the SC phase. 



The focus of the paper is the rather unusual properties 
of the SC phase. We will, therefore, assume throughout 
most of the paper that the SC phase does exist. This ap- 
proach is justified because it can under appropriate con- 
ditions be an equilibrium phase as discussed in Ref. jl^] 
and because even if it is not thermodynamically stable, 
there is a range of length scales over which correlations 
functions will exhibit SC behavior. 

This paper is composed of six sections, of which this 
is the first, and seven appendices, which mostly present 
mathematical details. Section O derives the Hamilto- 
nians for the columnar, sliding columnar, and nematic 
lamellar phases, based on a model in which the domi- 
nant interactions are couplings between each DNA lattice 
and the two membranes on either side of it. Section III 
presents the correlation functions for the SC phase, in- 
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eluding x-ray scattering intensities. Section [V addresses 



the stability of the SC phase. It calculates Td, and it 
derives general expressions for the interaction between 
dislocations before calculating Tkt- Section ^ presents 
a discussion of the various important length scales. Fi- 
nally, Sec. provides an overview of results. Appendix 
^ presents details of the derivation of the Hamiltonians 
in Sec. |l[ Appendices ^ through ^ provide details of the 
calculations of various SC phase correlations functions, 
and App. |^ presents details of calculations of interac- 
tions between dislocations. 



II. MODEL FOR LAMELLAR PHASES OF 
DNA-LIPID COMPLEXES 

A. Definition of Variables 

As depicted in Fig. |l|, lamellar phases of DNA-lipid 
complexes consist of a periodic stack of planar lipid bi- 
layer membranes with spacing a separated by galleries in- 
tercalated with DNA strands that form a local 2D smec- 
tic lattice with preferred spacing d. Since our primary 
interest is in the nature of possible ordering of the DNA 
strands, we will assume that the lipid membranes are free 
of defects such as dislocations or focal conic structures 
that destroy their integrity. We can, therefore, specify 
a layer by its integer layer number n, and we can spec- 
ify positions in the xz plane by the vector r = (x, z). 
We take the equilibrium height of the midpoint of bi- 
layer membrane n be no. The height of membrane n at 
position r is then 



(2.1) 



where /i"(r) is the Lagrangian height variable [|5[ mea- 
suring the displacement of layer n from its ideal height. 
The nth DNA lattice lies between membranes n and n-\-l. 
Its equilibrium height is, therefore, (n -I- i)a, and its 
height at position r is 



^^SnaW = {n+\)a + : 



:(r), 



(2.2) 



where (r) measure deviations from equilibrium height. 

The DNA lattice in layer n can be described in terms 
of a Fourier expansion of its density. It is important to 
keep track of the relative phases of mass-density waves in 
adjacent layers. In an ideal columnar lattice, the phase of 
the first mass-density wave in layer n is simply qo(-2 — z"), 
where 



qo ^ 2n/d 



(2.3) 



and z" specifies the preferred position of lattices in dif- 
ferent galleries. In a rectangular lattice, z" = 0, and in a 
centered rectangular lattice, z" = nd/2. Phase changes 
at position r relative to the above ideal phases are pro- 
duced by local translations described by the Eulerian dis- 
placement variable u"(r). Thus the DNA density- wave 
expansion in layer n is 



,ikqo[z~z"--u'^(r)] 



(2.4) 



If each DNA lattice is perfect with identical rods of linear 
density A and separation d, then p"'^(r) ^ X/d for every 
n and k. Thermal fluctuations will lead to reductions in 
p"'^ that are larger for larger k. Thus, p^NAi^) can be 
approximated by 



c.c. 



(2.5) 



where V', assumed to be independent of n, is the complex 
amplitude of the first non-trivial density wave. 

The total mass density of DNA strands at positions 

X = {x, y, z) = (r, z) is 

Pdna(x) = ^p"^-(r)e*'«"t^-^"-""("^^] 

nk 

x/^^a[2/-(«+5)«-<W], (2.6) 

where f^j^^ (y) is the form factor along the layer normal 
for a DNA mass-density wave of wavenumber kqo- If the 
DNA stands are lines with no width, the friNAiv) — ^iv)- 
If the strands are modeled as cylinders of radius ru^A, 
then the Fourier transform of f^^^A^v) 



foNAily) 



2Ji{^ql + k'^ql tdna) 



ql + k'^ql roNA 



(2.7) 



where Ji{x) is the Bessel function of order 1. The mem- 
brane density is 



Pmem(x) = ^ Pi /mem - UQ - /l"(r) 



(2.8) 



where p\ is the area density of the lipid and /mem(y) is 
the lipid membrane form factor, which typically is equal 
to 1/w for —1x1/2 < y < w /2 where w is the membrane 
width. X-ray scattering experiments probe (p(q)p(— q)), 
where p(q) is the Fourier transform of a linear combina- 
tion of pdna(x) and pmem(x) [See Sec. IIID|. 



B. Hamiltonian for Coupled DNA-Lipid Layers 

We are now in a position to derive the Hamiltonian de- 
scribing elastic fluctuations of lamellar DNA-lipid com- 
plexes. Our goal is to develop the simplest model con- 
sistent with symmetry. We begin with individual mem- 
branes and DNA layers. They are characterized by layer 
bending energies, which can be expressed as 

^K^a[{dlu-f + {dlK)% (2.9) 

where 0X3^ = Kmom is the bending rigidity of the lipid 
membranes and adK2d — kdna is the bending rigidity of 
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an individual DNA strand. We adopt a convention here 
with a factor of the layer spacing a multiplying sums over 
n to facilitate the continuum limit: ^ / '^^'^ ~^ I d^x. 
Interactions between DNA strands lead to a preferred 
separation between strands within a layer and to a com- 
pressional elastic energy, characterized by a modulus B2di 
for changing this separation. There is also a preferred 
separation between a given DNA layer and the two lipid 
membranes above and below it. Harmonic deviations 
from this separation are characterized by a modulus B^d- 
Changing the distance between membranes on either side 
of a DNA layer will lead to an expansion or a compres- 
sion of its DNA smectic lattice, which is described by a 
coupling with strength B^h between inplane strain and 
membrane separation. Combining all of these effects into 
a single compression energy, we obtain 



dxu 



n+l 



(2.15) 



where Ky « 4aV^''e"*<(**")''> . The interaction 7V^^ in Eq. 



( p'.13 ), favoring spatial registry, arises from interactions 
betw een DNA densities in Eq. (2^). The phases go-z" 
[Eq. (2.4)] are chosen so that energy is minimized when 
the remaining phase, /3"(r) = qo[z — u"(r)], of mass- 
density waves in neighboring DNA lattices are equal 
at the points of closest approach of the lattices [Fig. 
p. The point at r"+^ in layer n+l closest to the 
point at r" in layer n lies along the normal N"(r) « 
{~dxu];,l,-d,uf). Thus r"+i = r" (5r"+i, where 
(5r"+i = aN_L « {-adxu'^,0,-ad^u]^) where Nj. is the 
projection of N onto the xz plane. Energy is minimized 
when /3"+^(r"+^) = qo{z - adzUy - u""'"^) is equal to 
/?"(r"). The resulting energy is 



+ {Buh/a)u'^,{h 



n+l 



n+l 



where 



d.u:-mu:f + id.u:f]/2 



(2.10) 



(2.11) 



K - -V"" / d'^rcos[qo{u 



n+l 



is the nonlinear 2D Eulcrian strain and where factors of 



< + aa,0], (2.16) 



where is proportional to the square \tp\'^ of the mass- 
density- wave amplitude. As required, this energy is in- 
variant with respect to spatially uniform translations 
described by an n- and r-independent displacement of 
M"(r). It is also invariant to harmonic order with re- 



were introduced to facilitate the continuum limit, spect to uniform rotations in which u 



n+l ^ 



= -ad. 



e.g., (/i"+^ — h'^)/a ~ dyh{x.). There are additional com- 
pressional energies involving interactions proportional to 
(/i"+i - /i")2, (m^+1 - mJ;)2, _ u^j2^ and other 

further neighbor terms. These are smaller than those 
considered in Eq. ( |2.1Cl|) , and we will neglect them. They 
are, however, needed to achieve Td < Tkt and a stable 
SC phase [|3|. 

Neighboring DNA strands within a gallery prefer to be 
parallel. This leads to a Frank-like orientational energy 



The predominant interaction between DNA in neigh- 
boring strands is electrostatic, and we expect V"^ ~ 
d)e~'^^°-/'^ , where Ae is the charge per unit length 
on a DNA strand. 



1-L° 



(2.12) 



where ^"(r) « i9a;<(r) is the angle that the nth DNA 
lattice at r makes with the x-axis. In all but the NL 
phase, this orientational interaction is subdominant to 
those in Eq. (2.£), and we will ignore it. 

Finally, there are interactions between DNA lattices 
in neighboring galleries that favor parallel alignment and 
spatial registry of the lattices. These interactions are 
described by a sum of layer- coupling Hamiltonians 




The angular coupling is 



m = -V 



d2rcos[2(r-r+i)]. 



(2.13) 



(2.14) 



In all but the isotropic lamellar phase, there is long-range 
angular order, and we can replace by 



FIG. 3. A sequence of tilted and bent membranes. Interac- 
tions favor equality of the shifted phase, l3"{r) — qo[z — Un{r)], 
in neighboring layers at points of closest approach, i.e,, they 
favor /3"(r") = /3"+^(r"+i) where r"+^ - r" = aNx(r"), 
where Nx is the projection of the membrane layer normal 
onto the yz plane. 

To summarize, our complete model Hamiltonian cou- 
pling membrane height displacements /i" and DNA lat- 
tice displacements Uy and u" is thus 



n 



bend 



n 



int 



(2.17) 
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with entries defined by Eqs. (U), ( ^.lOD , and ( ^.13D to 



1. Effective Hamiltonian for DNA Displacements 

Though we will discuss membrane height fluctuations, 
our primary focus will be on properties of the DNA lat- 
tices. It is, therefore, useful to integrate out the mem- 
brane height variables h" (r) to obtain an effective Hamil- 
tonian for DNA displacements only. This operation is 
carried out in Appendix The long wavelength result 
is 



T/tot njSC 

r^DNA — 'T- 



(2.18) 



where W = ^„ HJJ and H.^'^ is the sliding phase Hamil- 
tonian, 

n •' 

-i^a j d\Ky[{d,u:+'~d,u:)/a]^ 



where 



r2 

J^zz — 



(2.19) 



(2.20) 



the summation convention on a = x, z and (3 = x, z is 
understood, and the bending-rigidity tensor is 



K 



a(3 _ 



K3d\ 



(2.21) 



Note that integrating out creates interactions between 
nearest-neighbor and next-nearest-neighbor u" and Uy 
displacements not present in the original model of Eq. 
(2.17). These interactions, however, have particular ra- 
tios determined by the parameters in 7i'°* . Of particular 
importance to us in what follows is the fixed relations 



between the coefficients of {Uy^ 



,n\2 



n+1 

zz 



<zY 



and u"^(uy+ — '). In a more general model with 
further-neighbor interactions, the simple relation of Eq. 
(2.20) among B^^-, B^h, and B^d would not hold. 

Using "H^Sna^ construct the long- wavelength 

Hamiltonian for each of the phases listed in the intro- 
duction as detailed in the following. 



2. The Columnar Phase 

The columnar crystal is characterized by strong cou- 
pling of displacements in different layers. Its elastic 
Hamiltonian can be obtained by expanding the cosine 
m W;, [Eq. (|I|)] about one of its minima. Performing 
this operation, taking the continuum limit, and retaining 
only the lowest order terms in a gradient expansion, we 
obtain the familiar elastic Hamiltonian for a rectangular 
columnar lattice |lql : 



BsdUyy + BuhUzzUyy + ByzUy^] 



+ I d^x[K2d{dlu,)^ + {K2d + K^dMuy)^] (2.22) 



where x = (r, na) 



(r), and B, 



V^q^a^. Note that the shear elastic modulus Byz arises 
from and is zero when is zero. The cross com- 
pression i3„/j-term arises from the 

in Ti}§\]A 1 which was generated by the interaction of it"^ 
with the membrane height field 



n+l 



ly) term 



3. Sliding Columnar Phase 

The sliding columnar phase is characterized by a van- 
ishing shear modulus Byz for relative displacement of 
DNA lattices. Its dominant fluctuations are, therefore, 
described by the Hamiltonian Ti.^'^ of Eq. ( 2.19| ) obtained 
by setting = in H^^^^. The coupling is irrel- 
evant when the sliding phase is stable as detailed in the 
following [See Sec. IV |. 



4-. Nematic Lamellar Phase 

Positional couplings between layers and the smectic 
compression modulus B2d both vanish in the nematic 
lamellar phase, and 



-fi [ cfx[B:iduly + +K''Hdc.do,Uyf]. (2.23) 



where = K2d and is a Frank elastic constant in- 
troduced in Eq. (|2.12|). 



III. 



CORRELATIONS IN THE SLIDING 
COLUMNAR PHASE 



In the preceding section, we derived a general Hamil- 
tonian capable of describing the phases of lamellar DNA- 
lipid complexes. In this section, we will investigate corre- 
lations in the sliding columnar phase. The most unusual 
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correlations in the sliding columnar phase are those in- 
volving displacements of DNA strands within the layers. 
We will thus begin by considering the effective Hamil- 
tonia n for u" obtained by integrating out Uy from H^^ 
[Eq. (2.1£)]. From this we will calculate all correlation 
functions of u". We will then consider membrane height 
correlations described by an effective Hamiltonian for ft," 
in w hich u" and u" have been integrated out of 7Y'°* 
[Eq. ( ^.17 )]. In the long- wavelength limit fluctuations in 



Uy are the same as those in ft". Finally, we will d iscuss 



the relevancy of the shear coupling V^, Eq. ( p. 16 ). This 
involves a consideration of both u" and dzU?, 



A. Correlations in it" 

1. Effective Hamiltonian for u" 

Since Uy and u" are harmonically coupled in H^j/^^, 
we can integrate over Uy exactly to obtain an effective 
Hamiltonian for u" fluctuations alone. This calculation 
is carried out in Appendix In the long-wavelength 
limit, the resulting Hamiltonian is 



where K = K2d and 



2d 



(BliJBsd)- 



(3.1a) 
(3.1b) 

(3.2) 



zz I 
SC 



Note that there is no term proportional to (u"+^ 
in Hamiltonian T^f*"^ even though there is one in H. 
This is a result of the special relation [Eq. ( ^.20 )] among 

energy coefficients in Ji^'-^ . 

The Hamiltonian of Eq. ( ^.l[ ) is the sum of two parts: a 
s um o f elastic energies for independent 2D smectics [Eq. 
( |3.1a)J and a term coupling angles in neighboring layers 
[Eq. ( 3.1b )]. The 2D smectic energy is a function of the 



full nonlinear strain w"^. When Ky — and there is no 
coupling between layers, nonlinearities in w"^ lead to im- 
portant renormalizations of the long-wavelength elastic 
constant of a 2D smectic [^ . When the interlayer cou- 
pling term is present, nonlinearities in vJl^ also lead to 
long wavelength rcnormalization of elastic constants in 
the SC phase [Q . These renormalizations are only loga- 
rithmic, however, and we will ignore them in this article. 

We, therefore, consider the harmonic limit of Tif^, 
which is conveniently expressed in Fourier space. Intro- 
ducing 



(2^)2 



^l'^ dqy dq^ dq^ 

-7r/a 27r 27r 27r 



where A^, and A^ ^ 27r/(i are wavenumber cutoffs, and 



<(r) 



d^q 



J(qyna+q±-r 



(27r)3 

where qj^ = (g^,, 0, gz), we obtain 



(3.4) 



^f''^^ / 4^[Bq'z+Kqt + Kyiqy)qlql]\uzi<ir. 



Here 



with 



(27r)3 



Kyiqy) = Kyp{qya) 



p{qya) 



^ (1 - cosgj^a) 
{qyaY 



(3.5) 



(3.6) 



(3.7) 



which tends to unity as qya ^ so that K{qy = 0) = Ky. 

Tw o important length scales can be obtained from Eq. 
( |3.l[ ) by comparing the orientational interaction energy 
with the 2D smectic compression and bending energies. 
The lengthscales 



— and z* — -^TT-, 
My My A 



(3.8) 



with /iy = y/Ky/K and A = y/K/B, separate two- 
dimensional from three-dimensional behavior. At length- 
scales within a gallery less than x* and z* , the 2D com- 
pression and bending energies are large compared to the 
orientational interaction, and the DNA lattices behave 
like independent 2D smectics. On the other hand, at 
lengthscales greater than x* and z*, the orientational in- 
teraction is significant, and 31? sliding behavior occurs. 

Before proceeding to discuss the sliding phase correla- 
tions, we n ote that the sliding phase elastic free energy 
in Eq. ( |3.l[ ) exhibits a striking local (gauge) translational 
invariance of the form 



Uz{x, z, na) Uz{x, z, no) + f{n)] 



(3.9) 



here the displacement /(n) is an arbitrary function of n 
assuming a different value in each gallery. In other words, 
the DNA lattices in different galleries can continuously 
slide relative to each other by arbitrary distances with 
no elastic energy costs . In a standard columnar phase, 
this continuous symmetry is spontaneously broken down 
by the shear coupling in Eq. ( 2.16| ). On the other side, in 
the sliding phase, this shear coupling is irrelevant at long 
scales, and elas tic properties reflect the gauge symme- 
try in Eq. (3.1). This symmetry is entirely responsible 



for the unconventional fluctuation behavior of the slid- 
ing phase we discuss in the following. In particular, we 
find that in the sliding phase, the fluctuations of relative 
displacements M""'"^(r) — u"(r) of DNA lattices in neigh- 
boring galleries diverge in thermodynamic limit. As a 
reflection of the gauge symmetry Eq. ([3.1|), a DNA lat- 
tice in a gallery is essentially free to flow relative to DNA 
lattices in neighboring galleries. The sliding phase thus 
exhibits zero macroscopic shear modulus [g|0] 
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S. Definitions of Correlation Functions 



Correlations in u" follow directly form Eq. (3.5) 



G..(r,na) = «(r)^.°(0)) ^ («(r))2) _ g(r, na) 



(3.10) 



where 



G..(q) 



Bg2 + + Kyqlqlp{qya) 



(27r)3 
and 

5(r,na) = i(K(r)-u°(0)]2) 



G..(q), 



(3.11) 
(3.12) 



d3 



G,,(q) 



1 _ g'i(qi-r+9y 



(3.13) 



It is useful to decompose g{r,na) into three parts: 

gir,na) ^ g^^^na) + g^^\r) - g^^\r,na), (3.14) 
where 

5W(na) = i([<(0)-^.°(0)n 
(fq 

j^:^Gzz{q){l - cos qyua) (3.15) 
5(^)(r) = i([u°(r)-zi«(0)]2) 

(2^ 

5(3)(r,n«) = -(K(0) - «°(0)]K(r) - <(0)]) 



/(g,G..(q)(l-e--) 



(3.16) 



x(l-cosgj^na)(l-e*i^''). (3.17) 

It is clear that g^-^'>{na) = g{0,na) and ^'■^•'(r) = g{r,0). 
Each of these functions has a characteristic singular be- 
havior as a function of system size and separation, which 
we will summarize below. Finally, the function 



g^'\r,na)^g^^\r)~g'^^\r,na) 



(3.18) 



will appear in our derivation in A pp. ^ of interplane den- 



sity correlations discussed in Sec. HID 



3. Asymptotic Forms for Ky 7^ 

The local fluctuation (('u")^) in the SC phase diverges 
as the square of the log of the system size with a func- 
tional form that depends on the order in which the sam- 
ple dimensions L^, Ly, and along the a;, y and z di- 
rections approach infinity: 



// n\2\ ^ 12 f lll^ [S>Lx/x*] Lz ^ Ly Lx 

U". ^ / - I 1 [a^Ljz*] » Ly ^ 
where the length 1^ is defined via 



T 



A^^^BKy 



(3.19) 



(3.20) 



ttz is a number, Lz ^ z* , ^ x* , and terms that do 
not diverge with system size have been dropped. The 
calculation of the displacement fluctuations in the limit 
Lz 00 and ~ Ly is detailed in Appendix Thus, 
SC "in-plane" fluctuations are less divergent than 2D 
smectic fluctuations that scale as a power-law with sys- 
tem size but more divergent than 3D Landau-Peierls 
smectic lamellar fluctuations that grow logarithmically 
with system size [p|JTo|]. The mean-square angular fluc- 
tuation, ((0")^) = {{dxUz)^), is finite, implying that the 
SC phase has three-dimensional long-range orientational 
order. 

The It? Lx^z divergence of ((u")^) is converted to a 
In^ r divergence in the function .g'^-'(r) [Eq. ( [3.16 )]. In 
Appendix y, we outline the calculation of g^^) (r) for large 
r. The results are 



.(2) 



(r) = ll 



2 / In^ [^e^x/x*] + Cx 



if z = 

i In^ [32e'^z/z*] + Cz if x = 0, 



(3.21) 



where 7 « 0.577 is Euler's constant and Cx and Cz 
are constants that depend on A^; and A^ but have well- 
defined Ax, 2 — + 00 limits. Note that the only Kx.z depen- 
dence of the correlation function in these two large dis- 
tance limits is found in the constants Cx and Cz ■ Cx and 
Cz are evaluated in Appendix ^ in the continuum limit 
{^x,z 00) in which we find Cx — and Cz — tt'^/S. 

The function g'^^^na) [Eq. (3.15)] is by definition zero 
when n = 0. For all n 7^ 0, it diverges logarithmically 
with system size: 



5W(na) = 2/2 5„(0) In 



An{Ax,Az 



where 



Sn{t) 



1 — cos(nit) 
du—= I, 
yjf^ + v?p{u) 



k=l 



2k -1 



(3.22) 



(3.23) 



(3.24) 



and An{Ax, Az) calculated in App. |^ depends on the 
layer separation n and the ultraviolet cutoffs but has a 
well-defined Ax^z ~> 00 limit that is calculated in Ap- 
pendix 

The function g^'^''{r, na) [Eq. (3.17)] is also zero at n = 
0. It has the same divergences as a function of separation 
r that (7^^) (no) has with system size: 
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5(3)(r,n)=2?X(0) 



HDnx/x*) 
ln(£'„az/z*) 



<lx 



-7/1/ 




0. 



(3.28) 



where -D„ and £"„ are numerical constants that have well 
defined values in the continuum limit, A^jA^ — > 00, as 
discussed in App. ^. 

It is useful to summarize the results of the calculations 
just presented. When n = 0, i.e., for points in the same 
layer, ^(r, 0) = (/'■^•'(r) grows with r as In^ r. When n ^ 0, 
i.e., for sites in different layers, g{r,na) diverges loga- 
rithmically with system size for all r because g'^^^ {no) di- 
verges in this way. If g^^^{na) is subtracted from ^(r, n), 
then the remaining function, g^'^\r,na), grows with r as 
In^ (r/r„) with r„ depending on the coefficient of Inr in 
g^'^'>{r,na) when n ^ 0. 



The scaling form of the correlation function implies that 
K2{qx,qz = 0) ~ ^nd ^2(9. = 0,g.) - 

In the harmonic regime where qx,zlx.z > 1, ^2(q±) — K2 
and i?2(q_L) = B2 are constants, and the scaling expo- 
nents are 7 = 4 and v = 2. In the anharmonic regime 
<lx.zix,z < Ij the scaling exponents 7 and u were cal- 
culated exactly by mapping the 2D smectic model with 
thermal fluctuations onto the KPZ model in 1 + 1 dimen- 
sions [To| ]. The exponents in the anharmonic regime are 
7 = 7/2 and i' = 3/2. 

The mean-square displacement fluctuations diverge in 
both regimes with lengths Lx and Lz of the sample in 
the xz plane: 



4- The limit Ky — 0; 2D Smectic Correlations 

When Ky = 0, there are no interactions between planes 
(when = 0), and the system reduces to a stack of 
independent 2D smectic planes. Since experiments ^] 
are consistent with nearly independent 2D smectic layers, 
we will review here well established results for such 
systems. 

Deco upled smectics are described by the Hamiltonian 
of Eq. (|3.1aD obtained by setting Ky = in H^^ . This 
Hamiltonian is a function of the full nonlinear 2D strain 
m"^. Nonlinearities lead to important deviations from 
harmonic behavior. Since there is no coupling between 
layers, Gzz{'^ is independent of qy, and we can define a 
2D correlation function (qj.) = Gzz(q)/a. 

At lengthscales less than the nonlinear lengths 



3/2 



T^/R2 



and /, 



A' 



(3.26) 



where K2 = Ka and B2 ~ Ba are 2D bend and com- 
pression moduli and A = \jK2lB2, nonlinearities are 
unimportant, and 2D smectic fluctuations are descr ibed 
by the linearized elastic Hamiltonian Tif of Eq. ( |3.5| ) 
with Ky = At lengthscales longer than and lz, 
the nonlin ear terms in the rotationally invariant strain in 
Eq. ( 2.1l| ) lead to renormalized bending and compression 
moduh K2{<\A_) and i?2(q_L) that, respectively, diverge 
and vanish at small wavenumber q_L Note that the 
nonlinear lengths lx and lz decrease with increasing tem- 
perature, and thus nonlinearities become important at 
high temperatures. In both the harmonic and nonlinear 
regimes, the Fourier transformed displacement correla- 
tion function (|M"(q)p) in each gallery were expressed 
as 



G2(q±) 



i?2(q±)g2+X2(q±)g^ 



where qx,z = qx,zlx,z and 



qi- 



(27r) 



G2(qi) = A2i2a/W(iz/L^), (3.29) 



where Uz{v) = M"(r), 2a = ^ — 1 — v = 1 in both regimes, 

Lx,z = Lx,z/lx,z, fi'\0) = const., and fi^\w) ^ 
as w ^ 00. This implies that the Debye- Waller factor 
(exp[i(7oUz])^ = exp[— gg(M^)] = in the limit of infinite 
system size, and there is no long-range positional order 
at any finite temperature in a 2D smectic, even when 
there are no dislocations. 

Since the mean-square displacement fluctuations di- 
verge as a power-law with system size, the displacement 
correlation function 



.2D 



(r) = i([<(r)-<(0)]2) 
^X'\x\'"fi'Hm/\m 



(3.30) 



diverg es algebraically with in-plane separation r. In Eq. 
( |3.30 ), X = x/lx, z = z/lzj and the scaling behavior of 

fu'\w) is similar to that of f^\w). In the harmonic 
limit, 



9o^5^^(r) 



where 



{z/^zf^ ifx = 0; 
{x/S.x), is z = 



XBld^ 
47r3T2 



(3.31) 



(3.32) 



is correlation length along z and — \{d / n)"^ {B2 / T) = 
2^/X£^Jtt is a correlation length along x. These correla- 
tion lengths are measured in x-ray scattering experiments 
at sufficiently short length scales when layers are decou- 
pled. 

Fluctuations in the angle 9 = dxU^ are nondivergent 
because of an additional factor of q^ in the numera- 
tor of ( [3.2£ ). Finite angular fluctuations imply that 



(3.27) {cosO) = exp[— (6'2)/2] is nonzero, and thus there is long- 
range orientational order in 2D smectics when there are 
no dislocations [pl. 
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B. Correlations in /i"(r) and Uy{r). 

In the absence of order in the DNA strands, lamel- 
lar DNA lipid complexes are simply multicomponent 
isotropic lamellar systems with height fluctuations iden- 
tical to those of any isotropic lamellar smectic. The pres- 
ence of orientational order in the DNA lattices introduces 
anisotropy into the effective bending moduli of the lipid 
membranes. By integrating out uj^(r) and M"(r) from the 
the complete Hamiltonian 7i*°* of Eq. (2.17), we obtain in 
App. ^ an effective Hamiltonian for height fluctuations. 
In the long- wavelength limit, this effective Hamiltonian 



Vh 



where 



=/(i^^7if^^if^7i^^^), (3.38) 



v^cos^ e + 2^i sin^ 9 cos^ O + v sin^ 9 ' 

(3.39) 



In Eq. ( 3.37 ), A and k{6) are cutoffs that depend on 
Aa;, Kz and ratios of bending moduli. In addition, A(6') 
depends on the angle 9 that r makes with the x axis. 



where 



B\am — Bsd — [BI^j-J B2d) 



(3.33) 



(3.34) 



and K"'^ is the bending-rigidity tensor defined in Eq. 
( ^.21 ). The model in Eq. (3.33) is an anisotropic ver- 



sion of the discrete smectic Hamiltonian often used to 
describe lamellar systems 18,1^]. An effective Hamilto- 
nian for DNA-lattice height displacements Uy{v) can be 
obtained by integrating over /i"(r) and u"(r) in 7i*°* of 
Eq. (2.17). The resulting Hamiltonian is identical in the 
long- wavelength limit to T-Lh in Eq. ( 3.33| ). Thus corre- 
lations in Uy{v) are identical in this model in the long- 
wavelength limit to those in h"{r). 

Fluctuations in h" are determined by the correlation 
function 



C. Density Correlations 

The DNA density correlation function arising from dis- 
placements parallel to lipid membranes is 



S{r,na) = (exp(igoK(r) - u°{0)]). 



(3.40) 



This function is easily evaluated in the SC phase when 
V = 0: 



S{r,na) = e-^«s(r,«a)^ 



(3.41) 



where ^(r, na) i s the displacement correlation function 
defined in Eq. (3.13). Since g{r,na) diverges with the 
system size for all n 7^ 0, S'(r, na) vanishes in the ther- 
modymic limit for all n 7^ 0. Thus, DNA densities in dif- 
ferent layers are completely uncorrelated in the SC phase 
when the coupling V" in Eq. ( 2.16| ) is set to zero. 



1. In-plane Correlations 



Ghh{q) 



(3.35) 



where K^f^qlqj = K^iql + qlf + K2dqt- This equation 
implies that ([/i"(r)]^) diverges logarithmically with sys- 
tem size as it does in ordinary isotropic lamellar systems. 
The coefficient of the system-size logarithm depends on 
the anisotropy in K"^ . Similarly, correlations in ft,"(r) 
defined via 

gn{v,na) ^ \{[h-(v) - h\0)f) 
= i([<(r)-<^(0)]2) 
d^q 



:Ghh{ci)[l 



_ g'(q-L-r+g!, 



(27r)3 

diverge logarithmically with na and r: 

Tjh \n{jiaXx^^^) if r = 



k^ghiv, na) 



2r]h ln(A(0)r) if n = 



(3.36) 



(3.37) 



When n = 0, g{r, na) does not diverge in the thermo- 
dynamic limit, and we have 

5*2(0 = Sir,0) = (e^'oKW-^'Wl) = g-'^f'^W, (3.42) 



with o^r) given by Eq. (|1|). Thus, from Eqs. ( |3.21| ) 

and 



52(x,0) = ^-y. ' ^>>^; (3.43) 

L e 'off i^^uj^ a; < x*, 

where Sx = e~'^o^^'^^ ig a constant, and g^^{r) is the 
displacement correlation function of a 2D smectic. In 
the other direction, 

S2i0,z) = \^^' ,^ ' ^>>^; (3.44) 

where — e^'^o^"^' is a constant. 



where A^, = ^/B/K^, fcp = 27r/a, and 
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2. Inter-plane Correlations 

As we have just seen S{r,na) is zero when n ^ if 
= 0. When is not zero, S{r, na) is not zero, 
and it can be calculated perturbatively in a ex pansion in 
V^/T. Using the decomposition of Eq. ( pl^ ) of Hg^r^ 
in to Ti^'-^ and 7i", we can express S{y, na) as 



5(r, Tia) 



^p--H"/TgZ?oK(r)-«°(0)]^ 



(3.45) 



where (•) signifies an average with respect to H.^'-^ . This 
expression is easily expanded in a power series in /T. 
The evaluation of even the lowest order term in this ex- 
pansion cannot be carried out analytically. In App. ^, 
we evaluate the lowest order term in an approximation 
in which Ky = 0, But = 0, and correlations in dzUy{r) 
in different layers are ignored. The result for the Fourier 
transform of S{r, na) is 



S{q±,na) 



^) [^2(q±)]"+\ 



where = V^'e'^y with 



Wy=qy{{d.u-y)')/2 



and 



S2iq±) 



(27r)2 



""52 (r) 



(3.46) 



(3.47) 



(3.48) 



is the Fourier transform of the two-dimensional intra- 



plane density correlat ions function of Eq. (3.42). Since 



the expression, Eq. (3.46), for S{q±,na) was derived 
under the assumption that Ky = 0, S'2(r) is strictly 
speaking the density correlation of an isolated 2D smec- 
tic. However, a more sophisticated variational approx- 
imation, to be presented in a separ ate p ublication [ p^ , 
yields a result very similar to Eq. ( 3.46 ) for S{ci±,na) 
with S'2(q_L) replaced by the true intra-plane correlation 
function with Ky ^ obtained by Fourier transforming 
5(r,0) of Eq. ( |3^ ). 

Transforming S{q±, na) back to real space, we obtain 



Sir,na)^(pi J 



(2^)2 



g2q-L-rg(n+l) In S2{q±} 



(3.49) 



The large n limit of this expression can be calculated 
using steepest descents. For fixed r and n — > oo, the 
result is 



S'(r, na) 



2nT 



where 



Sa 



S2{r = 0) 



d\rlS2{r) 



(3.51) 



for a = X, z, and Tq 



^j, = a In 



X, z and 
2T 



V'-S2iq± = 0) 



(3.52) 



Since 5*2 (r) dies off at least as fast as e when Ky = 
and as e xp(— const. In^ r) when Ky ^ 0, the 



tegrals in Eq. ( ^3.51 ) converge, and £,x and ^2 are well 
defined. 



D. X-ray Scattering 

X-ray scat- 
tering experiments probe /(q) = (p(q)p(— q))/^, where 
V is the volume of the system, and p(q) is the weighted 
sum of the DNA and membrane densities poNAi'^) and 
Pmcm(x). The total scattering intensity can thus be bro- 
ken up into a DNA contribution loNAiq), a membrane 
contribution /mcm(q), and a DNA-membrane cross term. 
Scattering from the ordinary columnar phase will exhibit 
true Bragg peaks at G;™ — (0, {I + ^am)kQ , mq^) . Here 
(7 = for a simple rectangular columnar lattice, and 
a — \ for a centered rectangular columnar lattice. Since 
there is no long-range positional order in the SC phase, 
there will be no Bragg peaks in its x-ray scattering profile. 
Rather there will be membrane-dominated lamellar peaks 
at G/o = (0, Zfco, 0) and DNA-dominated sub-power-law 
peaks at Gim for m 7^ 0. 

The contribution of lipid membranes to the x-ray scat- 
tering intensity is 

1, 



-^incm(Q) — I /mem (O'y ) I 



X Ve-*«»'"^ / (fre-'''^-'Sh{r,na,qy), (3.53) 



where 



5„(r,na,g,) = (e-'«='['*"«-''°(°)l) 



— g-9„ffh(r,na) 



(3.54) 



and 5h(r, na) is evaluated in Eq. (3.36). Thus, for exam- 
ple. 



Sh{Q,na,ko) ^ (na)-'"^ 



(3.55) 



Except for some anisotropics, this yields power-law Bragg 
peaks that are essentially identical to those of a normal 
lamellar smectic. 

The contribution of the DNA lattices to the x-ray scat- 
tering intensity is 



lDNA[q) = {PDNA{q)PDNA{-q)) 



-ikqo 2" 



nk 



d^re-'^i' •'"S'y2(r, na, qy, k), (3.56) 
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where Aq^^ = q± ^ kgoe^ and 

Sy,ir,na,qy,k) = (e-*««KM-<(o)]e-»9ofcK(r)-«°(o)]^^ 

(3.57) 

This correlation function is even more compUcatcd to 
calculate than S{r,na) [Eq. j^A^]. When B^h = 0, u"^ 
and decouple, and when fc = ±1 we have 

Syz{r,na,qy,k = ±1) = Sh{r,na, qy)S{r,na), (3.58) 

where we used the fact that correlations in Uy are the 
same as those in /i". Thus, at r = 0, this function dies 
off as n~'"'~^e~"°/'5!' when qy — fco- 

If we ignore the Sh{Y, na, qy) contribution of Syz (which 
would be justifie d for i?iam — > oo) and we use the approx- 
imation o£Eq. ( [3.49| ) for S{r,na), then the sum over n 
in Eq. (3.56) can be carried out exactly. The result for 
the dominant contribution with /c = ±1 is 

lDNA{ci) = i^2(Aql)|/A^^(g,)|2F(q), (3.59) 



where 



1 - 



2T 



1 _ S ^-^-f'^'^ cos(g,a - a.) + (^^^^^ 



(3.60) 

The structure function F{q) reaches a maximum at 
qy = mfco if a rectangular (cr = 0) lattice is preferred 
and at qy — {m + ^)ko when a centered rectangular lat- 
tices (cr = 1) is preferred. At these peaks, F{q) does 
not exhibit singular behavior: the function S'2(Aqj^) [Eq. 
( p8| )] entering Eq. ( ^!60| ) can be expanded in powers of 
Aqj^ to all orders. 

Ou t-of-pl ane lamellar fluctuations, ignored in Eqs. 
(|^ and ( |l60|) (by assuming Sam ^ oo), become sig- 
nificant when the inter-layer positional coupling V" is 
weak. For exam ple, for = and finite i?iam, we find 
using Eqs. {KE^ to (^!58|), 



IdNA - - \fDNAilv)\'^ 

X J d'^re-'^'^'^ ''Sh{r,na^0,qy)S2{r), (3.61) 



Using Eq. ( p. 541 ), it is easy to see that Idna{(i) has a 
maximum at qy = 0, Aq^ = ( even for infinitely thin 
DNA, with f]jNA{qy) — 1). This peak remains at gj, = 
for sufficiently weak nonzero jlij. Thus, for suffi- 
ciently weak interlayer coupling, the scattering structure 
factor has a form resembling that of a simple rectangular 
lattice, with a (0, l)-like peak at gj, = even though the 
system has short-range centered rectangular order in real 
space, With increasing strength of inter-layer positional 
coupling, this peak will bifurcate, at a critical value of 



V", into two (±1,1) peaks at nonzero qy [as in Fig. |j. 
Such a change of the form factor is not accompanied by 
a thermodynamic phase transition. It is similar in char- 
acter to the transition ac ross so -called disordering line in 
random microemulsions 



21 1 at which the wave- vector 



maximizing S'(q) vanishes. 



IV. STABILITY OF THE SLIDING COLUMNAR 
PHASE 

Until now, we have treated the sliding columnar phase 
as though it were thermodynamically stable. We will 
now examine conditions for this stability. In order for the 
SC phase to exist, it must be stable against forces that 
bring about registry between neighboring DNA lattices 
to produce the columnar phase; and it must be stable 
with respect to the disordering effect of dislocations that 
leads to a nematic lamellar phase. In this section, we 
will investigate these two effects. We find that the SC 
phase orders into the columnar phase via a roughening- 
like transition at a decoupling temperature Td and that 
it melts to the nematic phase via a dislocation unbinding 
transition at a temperature Tkt- Thus, the SC phase 
can only exist in a temperature range Td < T < Tkt, 
and it cannot exist at all if Tkt < Td- We calculate 
Td and Tkt for our nearest-neighbor model, and we find 
that 



Td TT^ 



(4.1) 



for systems with temperature independent coupling con- 
stants. Thus Tkt < Td, and, for our model, the SC phase 
is never thermodynamically stable. The introduction 
of competing nearest-neighbor and next-nearest-neighbor 
strain couplings can, however, reverse the inequality on 
Tkt and Td and stabilize the SC phase jl^]. Even when 
the SC phase is not thermodynamically stable, there is 
a range of length scales in the lamellar nematic phase in 
which correlation functions will exhibit SC behavior. 



A. Relevance of Translational Coupling 

In the sliding columnar phase, fluctuations of relative 
displacements M""'"^(r) — u"(r) of DNA lattices in neigh- 
boring galleries grow with increasing systems size. We 
found in S ec. |IIIA2| that ([<+Hr) - ""(r)]^) - InL 
[See Eqs. ( 3.15 ) and ( |3.22 )]. Due t o this divergence, 
the translational coupling [Eq. ( ^.16 )] may become ir- 
relevant above a critical decoupling temperature. This 
result is obtained by calculating the expectation value of 
the translational coupling 



{ni) = -V-J <fr{cos[qoii 



n+l 



u- + adzum (4.2) 
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with respect to the shding columnar Hamiltonian in ( 3.1 ). 
Because the cross correlation {dzUy{u^'^^ — u")) is zero, 
we can evaluate this quantity exactly: 



-w. 



J d^rexp[- 



(4.3) 



where Wy = q^a {{dzUyY) /2 , Using the fact that 
i) diverges logarithmically with system size [Eq. 



( 3.22 )], we find that the translational coupling scales as 



(K> - -y^e-^^i^-r,.^ 



where 



AT 



(4.4) 



(4.5) 



The condition 2 
temperature 



lid 



<P^BKy 
defines the critical decoupling 



d^^BKy 



(4.6) 



When T < Td, the translational coupling scales as sys- 
tem size to a positive power, and is relevant. In this 
case, the system becomes a columnar phase at the longest 
lengthscales with a nonzero shear modulus for shifting 
neighboring lattices relative to each other and long-range 
positional order in the yz plane. When T > Td, the trans- 
lational coupling scales as the system size to a negative 
power, and is irrelevant. The system flows to the 
sliding columnar phase at the longest lengthscales when 
T > Td- Thus, Td marks the transition from the colum- 
nar phase to the SC phase as temperature is increased. 
This transition is of the roughening type. There is no 
energy cost for shifting neighboring lattices relative to 
each other, and thus the sliding columnar phase is po- 
sitionally disordered in the yz plane. These conclusions 
are supported by an RG analysis to be presented in Ref. 



B. Dislocation unbinding 

We have just seen how thermal fluctuations weaken 
interlayer couplings to produce the SC phase from the 
columnar phase. We will now consider the disordering 
effects of dislocations. An individual edge dislocation in 
a 2D smectic has a finite rather than a logarithmically di- 
vergent energy. As a result, there are thermally excited 
unbound vortex pairs at all temperatures that convert 
the 2D smectic into a 2D nematic at the longest length- 
scales [1). In a sliding columnar phase, each DNA smectic 
layer experiences an orientational ordering field from its 
neighbors. As a result, the energy of an individual edge 
dislocation in a given layer diverges logarithmically with 
system size, and there can be a Kosterlitz-Thouless dislo- 
cation unbinding transition. The low-temperature phase 
with bound dislocations is the sliding columnar phase, 
and the high-temperature phase with unbound disloca- 
tions is the lamellar nematic phase. 



1. Dislocation Energy 

The DNA smectic lattice in each layer can have edge 
dislocation defects in which the displacement field un- 
dergoes a change of kd^ where k is an integer, in one cir- 
cuit around the defect core. If there are dislocations with 
integer strengths fc„^; at positions r„^; in layer n, then 



X v"(r) = 6"(r)y, 
where v"(r) = Vj^u"(r) and 

V'{v)^dY,kr,d5\v-Vr,^l) 



(4.7) 



(4.8) 



is the dislocation density in layer n. 

We can now calculate the energy cost of an arbitrary 
assembly of edge dislocations in the sliding columnar 
phase using the SC Hamiltonian 



d^r 



B[v:f + K{d.vlf 



2 ^ X 



,,n+l\2 



(4.9) 



written in terms of v"(r) = (i^",'*^")- Our strategy is to 
minimize this Hamiltonian subject to a nonzero disloca- 
tion density 6"(r). As usual, the calculation is simpler in 
Fourier space. To transform from real space to Fourier 
space, we use 



dVe-^^i^ "'')v"(r), 
6(q)=^a J d2re-*('i^""+«''™)6"(r) 



ad^^ kn. 



le 



(4.10) 



(4.11) 



71,1 



We minimize Eq. (^^) using the Euler-Lagrange equa- 
tions and flnd 



B 



where 



K{c)q^ ^Kql+Kyqlp{q, 



(4.12) 



(4.13) 



We then employ the constraint, Eq. ( |4.7| ), to relate Vz(q) 
to the specified dislocation density fo(q). In the final step, 
we insert the expressio ns fo r v(q) into the Fourier trans- 
formed version of Eq. (4.9) and find 



Ed — 77 



77 1 a 



7r/a 



dq^ f d^q^ K{q)q^\b{q)\^ 
27: J (27r)2 Sg2 + Kic^q^ql 



(4.14) 



for the energy of edge dislocations in the sli ding c olumnar 
phase. Also note that if we set Ky = 0, Eq. ( 4.14 ) reduces 
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to the expression for the energy cost for edge disfocations 
in a 2D smectic §. 

We now decompose the disfocation energy into a part 
that diverges with system size and a part that diverges 
only with the separa tion b etween defects. After we insert 
Eq. ([4. Ill) into Eq. ([4.14|) , we find 



= ^^nf^n' Jn-7i' In [B^-n' i^x) / X*] 

n.n' 



.1,1' 



(4.15) 



where tT„ = 
7ith layer, 

Jn, — 



kn,i is the total dislocation charge in the 



du^/v?p{u) cos nu — — 4~2' (4.16) 



and Bn{^x) is given in Appendix |^. The first term di- 
verges logarithmic ally w ith system size. The interaction 
energy £'„(r) [Eq. ( [1.15| )], on the other hand, is written in 
terms of an integral that vanishes when r = 0, and thus it 
diverges only with the separation between defects. In the 
limit of infinite system size, Jn-n' is a positive definite 
matrix. The dislocation energy diverges logarithmically 
with system size unless the total dislocation charge in 
each layer is zero, i.e., cr„ = for every n. 



The dislocation interaction energy of Eq. (4.15) 



dt 



£;„(r) = - du —\/t^ + u^p{u) cos nu (4.17) 

Jq Jo i 

X ^1 — cos[tx/x*] exp + u^p{u) ^ , 

is difficult to calculate for arbitrary separations (r, na), 
however, it can be calculated in the limits x ^ x* , z = 
and z :s> z* , X = 0. We find that £'„(r) scales as 



E„ 



r ln[C-(A,)|x|/x* 
-•'-\HC^{A,)\z\/z* 



if X 3> X* and z — 
if z ^ z* and a; = 0, 



where C!^''^{Ax) are numbers that depend on the layer 
index n and the cutoff A^; and are given in Appendix |G|. 
Note that the A^, ^ oo limit of £'„(r) is not well-defined. 
Since £'„(r) ~ lnr/(4n^ — 1) for large r, -E„(r) is positive 
for all n > and negative for n = 0. As a result, like- 
signed dislocations in different layers attract each other, 
whereas like-signed dislocations in the same layer repei 
each other. 

The attraction of like-sign dislocations in different gal- 
leries makes physical sense. The dislocation excitations 
of the SC phase can be viewed as closed loops carrying 
a single value of charge with portions of the loop passing 
normal to layers and portions passing parallel to layers. 
Those parts of the loops passing through layers give rise 
to layer dislocations, and those parts aligned parallel to 
layers cost no energy because the shear modulus is zero 



[^,0. A direction can be assigned to a loop. A loop sec- 
tion penetrating a layer in the upward direction gives rise 
to a dislocation of one sign while one penetrating in the 
downward direction gives rise to a dislocation of the op- 
posite sign. As is apparent from Fig. ^, the expectation is 
that the lowest energy configuration of a continuous loop 
will penetrate nearby points in nearby layers in the same 
direction, i.e., that nearby dislocations in nearby layers 
will have the same sign. This argument is, of course, 
not rigorous because it is possible for dislocation lines 
running parallel to layers to cross to produce nearby dis- 
locations in nearby layers of opposite sign as shown in 
Fig. ^(e). In fact, it is possible to construct interlayer in- 
teractions for which like sign dislocations in neighboring 
layers repel rather than attract. 



+ 






+ 






> 

+ 







(a) 




2+ 2- 



2±. 



(c) 



(d) 




(e) 

FIG. 4. (a) A sequence of closely bound dislocations in 
neighboring layers interpreted as a dislocation loop, (b) An 
unbound version of the dislocation configuration of (a), also 
interpreted as a dislocation loop, (c) a complex compos- 
ite-dislocation pair, (d) Melted composite pair of (c). (e) 
Opposite sign dislocations at nearby points in neighboring 
layers interpreted as a crossed dislocation loop. 



2. Dislocation Unbinding Temperature 

In the previous section, we showed that the dislocation 
energy is logarithmically divergent unless there is dislo- 
cation charge neutrality in each layer. Consider now a 
configuration in which layer n has a single dislocation of 
charge a„ (which may be zero). The free energy of su ch 
a composite dislocation configuration is, by Eq. (|4.15| ), 



Fn = En- TS 



D 



(4.18) 
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Jn—n' 



2T lnL + const. 



Clearly Eq. ( 4.1^ ) indicates that composite dislocations 
pairs unbind for temperatures above the Kosterlitz- 
Thouless type temperature, 



TKT[<7n] ~ 



87r2 



n.n' 



(4.19) 



If cr„ is nonzero only in layer m (cr„ — Snm), the transi- 
tion corresponds to unbinding of simple dislocations pairs 
within a single layer. If cr„ is nonzero in more than one 
layer, the transition corresponds to the unbinding of com- 
posite multilayer dislocations pairs as depicted in Fig. 

We see that each configuration of dislocations {<Tn} 
will have a different unbinding temperature. The naive 
expectation is that Tkt is lowest when there is a sin- 
gle dislocation in a single layer. It is, however, possi- 
ble that composite dislocations might melt at a lower 
temperature. Hence, we define Tkt — niinCT„ rfj-T[(7„]. 
Above this temperature, one of the dislocation configu- 
rations (either composite or individual) unbinds, renor- 
malizes the compression modulus B to zero, and destroys 
the in-plane 2D smectic order of the sliding columnar 
phase. The lowest unbinding temperature corresponds 
to the configura tion with t he lo west Ini^ divergent en- 
ergy [See Eqs. (|]l|) and (p^ j. Since J„ < for aU 
n > 0, \Jn\ decays with increasing layer separation, and 
the dislocation energy scales quadratically with the dis- 
location charge, the lowest energy configuration with N 
defects is a string of N defects of strength -1-1 with one 
in each layer, i.e. (T„ — ^n,p- It is then straightfor- 

ward to show that the lowest of these configurations is an 
individual defect with = 1. Therefore, the dislocation 
unbinding temperature for the sliding columnar phase is 



Tkt — 



27r2 



(4.20) 



The ratio of t his t emperature t o th e unbinding temper- 
ature T^ [Eq. (0)] yields Eq. {1^ for (3 = TKT/Td = 



C. Comments on Lyotropic Systems 

An interesting aspect of the present study is the tran- 
sitions from the sliding to the columnar and the nematic 
phase we discussed above in terms of the correspond- 
ing critical temperatures: the decoupling temperature Td 
[Eq. ( |4.6| ], for t he co lumnar to sliding phase transition, 
and Tkt [Eq. ( 4.20 )] for the sliding to nematic phase 
transition. In lyotropic systems of interest here [Refs. 
[^- ||], temperature changes can hardly produce signif- 
icant effects. Nonetheless, the above transitions can be 
still triggered by varying the period of DNA lattices (as 



lattices described in Refs. [g|). Indeed, by previous 
discussions, the decoupling condition T > Td, is equiv- 
alent to rid > 2, with r]d = AT / (P{BKyY/'^ . Thus, the 
decoupling transition from the columnar to sliding phase 
may be triggered by changing the strength of the inter- 
layer orientational coupling constant Ky. As evidenced 
clearly in experiments the strength of inter-layer 

couplings significantly increases by increasing the period 
of DNA lattices. Thus, by the swelling these lattices, 
one may reach the transition in which the sliding phase 
changes into a columnar phase, as discussed in Sec 



[VA 



and, in more detail, in Ref. [ |14| . Finally, we emphasize 
that, within the simple model discussed here, there is 
no true shding phase, as suggested by Eq. (11) that ap- 
plies to systems with temperature independent coupling 
constants. It should be stressed however that the restric- 
tion of having temperature-independent constants is not 
essential for this important conclusion about the slid- 
ing phase stability. Neither it is essential whether the 
system is thermotropic or lyotropic in its character. In- 
deed, by Sec. tV A , the condition to have dislocations 
bound [i.e., T < Tkt] is equivalent to ijkt > 2, with 
rjKT = (P{BKyY/'^ /tt'^T , for the model discussed here 
[Sec. 0. Note that 



Vd Vkt = A/tt^ < 4. 



(4.21) 



This relation implies that the condition for dislocations 
confinement {tikt > 2) and the condition {rjd > 2) for 
sliding phase decoupling cannot be simultaneously re- 
alized. Whenever the decoupling condition is realized, 
rjd > 2, rjKT = A/n'^rjd < 2, and dislocations will unbind. 
This feature turns the sliding phase into a nematic phase 
at long length scales. Still, at low dislocation densities 
(what may well be the realistic case, see Sec. in such 
a nematic phase, there is a broad range of length scales 
exhibiting sliding-phase correlations and other structural 
properties discussed in this work. 



V. CHARACTERISTIC LENGTHSCALES 

Whether or not sliding columnar behavior is seen in re- 
cent X-ray scattering experiments on CL-DNA complexes 
is still an open question Even though the slid- 

ing columnar phase is converted into a nematic lamellar 
phase at the longest lengthscales, the SC phase may exist 
at shorter lengthscales determined by the density of edge 
dislocations. Thus, CL-DNA complexes studied in recent 
experiments with domain sizes L « 0.1//m may exhibit 
sliding columnar behavior. However, in the next round 
of experiments it will be important to prepare aligned 
CL-DNA samples because powder-averaging complicates 
the functional form of the scattering intensity and makes 
it difficult to identify sliding columnar behavior. 

We showed in Sec. Ill that the density-density corre- 
lation function 5*2 (r) [Eq. ( |3.42| )] displays different func- 
tional forms depending on the magnitude of the in-plane 
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separation r. The crossover lengthscales for the correla- 
tion function are l^, Iz, Cd, x* , and z* . The harmonic 2D 
smectic regime is defined by x, z < l^^z and the nonhnear 
2D smectic regime is defined by x, z > l^.z, where the 
nonhnear lengths l^^z are given in Eq. (3.26). To esti- 
mate the nonlinear lengths, we must determine and Iz 
in terms the experimentally measured qu antiti es d and 
the harmonic correlation length of Eq. ( 3.32 ). and 
ly, can be written in terms of d and ^z by solving Eq. 
32) for i?2 and using 



(1 



K2 = np/2d, 



(5.1) 



where = 500A is the persistence length of DNA. Ta- 
ble I evaluates the nonlinear lengths for two DNA spac- 
ings of the experiments of references [^J^ . It shows that 
L < Ix.z and indicates that significant departure from 
harmonic 2D smectic behavior is not expected in agree- 
ment with these experiments. 

The finite length of the DNA molecules ^na ~ 16/im 
introduces another crossover lengthscale. The density of 
DNA molecules within a given layer is p = 1/ d^oNA- If we 
assume that each free end of a DNA strand corresponds 
to a dislocation, then we can estimate the characteristic 
distance between dislocations to be 



DNA- 



(5.2) 



This length is clearly smaller than the actual distance 
between dislocations because not every end of a DNA 
strand has to produce a dislocation: a series of DNA 
strands can align in a row to produce a layer of a 2D 
smectic. At lengths scales greater than ^d, the 2D 
smectic behaves like a nematic |^. Using the estimate 
of Eq. (5.2), we find « 0.21/zm and 0.30/im when 
d ~ 28A and 55 A, respectively. Note that L < S^d, 
and thus the subdomains are small enough to possess 
2D smectic ordering. 

TABLE I. The nonlinear lengths, and Iz , calculated as a 
function of the DNA spacing d using the experimental values 
of the in-plane correlation length . 



Nonlinear length 


d = 28A 


d = 55A 




0.35/im 


0.24^im 


h 


9/im 


6/im 



We can also estimate the energy cost for creating 
hairpin edge dislocations within the 2D smectic lattices. 
Hairpins cause the DNA director to change by tt over a 
lattice spacing d. The energy cost for a hairpin can be 
estimated from the 2D bending energy. If we ignore relax- 
ation of the 2D smectic layers, we estimate the bending 
energy to be 



Ehp 
2T 



4 d' 



(5.3) 



which implies that hairpins are favored on lengthscales 
greater than 



ihp = d exp 



'2T 



(5.4) 



Since (^hp ^ > L throughout the experimental range 
in d, hairpins are not important for the current set of 
experiments. We do not yet have accurate estimates of 
X* and z* since the value of the orientational rigidity 
Ky is unknown. Scattering experiments will see sliding 
columnar behavior on lengthscales less than £^d if x* ,z* < 
id. 



VI. DISCUSSION AND CONCLUSION 

In this paper, we have introduced the new sliding 
columnar (SC) phase of matter which may exist in lay- 
ered systems composed of weakly-coupled 2D smectic 
lattices. The sliding columnar phase is characterized 
by weak positional but strong orientational correlations 
between neighboring 2D smectic lattices. The SC har- 
monic free energy contains an orientational rigidity that 
aligns neighboring 2D smectic lattices in addition to in- 
plane compression and bending moduli. The SC phase 
is characterized by a vanishing shear modulus for rela- 
tive displacements of 2D smectic lattices. The presence 
of the orientational rigidity fundamentally alters the en- 
ergy spectrum. In light of this, we have calculated the 
structural properties of the sliding columnar phase such 
as the SC displacement correlation function, scattering 
intensity, and dislocation energy. 

Experimental research on layered liquid crystalline sys- 
tems studied in this work is in progress §- §. Corre- 
lation peaks theoretically discussed here have been ob- 
served in a number of scattering experiments on DNA- 
cationic lipid complexes. Reference reports DNA 
scattering patterns that clearly reflect the short-range 
centered-rectangular order depicted in Fig. |[ Reference 
1^, on the other hand, reports peaks in a different sys- 
tem at the simple rectangular lattice positions (0, 0, ±qo) 
[i.e., at (0, 1) peaks]. In spite of this difference in scatter- 
ing data, it is likely that both systems have short-range 
centered rectangular order in real space. Indeed, in a 
strongly fluctuating disordered systems, positions of cor- 
relation peaks in g-space may not always reflect length 
scales and short-range order in real space. For example. 
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the correlation peak in g-space in random microemulsions 
po| , pl| may occur at g = even in situations in which 
real space data still show a fin ite st ructural length scale. 
As discussed at the end of Sec. [II C| , for sufficiently small 
coupling inter-layer positional coupling V""/r, the scat- 
tering form factor loNAi^i) exhibits a maximum at the 
simple rectangular-lattice points (0,0, ±(jo) even though 
the system has short-range centered rectangular order in 
real space. The differences between the positions of the 
scattering peaks in Refs. |^ and may be explained 
by the differences in thermal disorder and/or inter-layer 
couplings strengths of the two systems. The interlayer 



is several times 
The former sys- 



correlation length reported in Ref. 
longer than that reported in Ref. 
tem is stiffer and more ordered than the latter. Conse- 
uently the (1, 1) and (—1, 1) correlation peaks in Ref. 
^ reflect the centered-rectangular order in real space 
depicted in Fig. |l|, whereas the the (0, 1) peak reported 
in Ref. |^ reflects the merging of (±1,1) peaks brought 
about by strong thermal fluctuations and/or sufficiently 
weak inter-layer couplings. 

After completing this work, we learned from Joachim 
Radler of experiments p2| showing a continuous varia- 
tion of the DNA form factor with changing strength of 
inter-layer coupling consistent with the change we 
describe above and in Sec. [IIP [crossover from centered 
rectangular scattering pattern to one resembling that of 
a simple rectangular lattice]. In these experiments, the 
variation of inter-layer positional coupling is produced by 
varying the period of DNA lattices in galleries. This cou- 
pling generally increases with increasing period of DNA 
lattices, most likely because the Coulomb interaction in- 
creases and the elastic membrane deformation increases, 
as suggested by electron density maps of the system [p2[ . 
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APPENDIX A: DERIVATION OF EFFECTIVE 
HAMILTONIANS 

In this append ix, w e will derive the effective Ha milto - 
nians, [Eq. (|2A|)], Hf^ [Eq. (P)], Hh [Eq. ( p3| )] 
obtained, respectively, by integrating out /i"(r), ft,"(r) 
and u;'(r), and <(r) and u^(r) from H = ^•^'="'^+7^™'"-!- 
7^"^°*, which we can express in linearized form in Fourier 
space as 



(27r)3 



- ^[A,(q)^i,(q)/i(-q) + A,(-q)u,(-q)/i(q)]}, (1.1) 

(T 

where a — y^z^ and 

A(q) =i^3dgi+4(S3d/a') 

Ay{ci)=A{BM/a')+K2dqt 

^^(q) = B2dql + K2dqt + Kyqlqlp{qya) 

\y{ci) = (2B3d/a2) (1 + e-«9.a) 

\M) = -{Buhla)iq, (e-'««'^-l). 

Integration over ft.(q) yields 



(1.2) 



use 



(2^ 



where 



-//(q)Mj^(q)M;j(-q) - ^(-q)uy(-q)M3(q)], (1.3) 
|A.(q)P 



M(q) = 



In the small qj^ limit, these expressions reduce to 
4(q) = (B3d/a')2(l - cosg.a) + K'-^qlql 

A;(q) = B2dql - ^qlni - cosg^a) + K2dqi 



3d 



+Kyqlqlp{qya) 
/i(q) = -{Buh/a)qzSini 



(1.4) 
(1.5) 

(1.6) 



(1.7) 
(1.8) 



where K"^^ is the bending rigidity tensor of Eq. (2.21). 
Fourier transformation of Eq. ( |l.3| ) back to real space 
using the low-qj^ form of Ay (q) , A^ (q) and /i (q) p roduces 
the sliding phase Hamiltonian TY'^'^ of Eq. ( |2.19[) . 

The effective Hamiltonian for is obtained by inte- 
grating 'H^'-^ 



over Uy. The result is 



1 



Sq 
(2^ 



where 



A';(q)=A;(q) 



lM(q)P 



(1.9) 



(1.10) 



In the q_L — > hmit, ^"(q) becomes 
A"M) ^ qt B 



2d 



-^2(1 



i%a) 



sm qytt 



B3d 2(1 - cos qy a) 
+Kyqlqlp{qya) + K2dqt 



(1.11) 



B2d — 



B 
B 



uh 



3d 



K2dqxi 



+Kyqlqlp{qya) 



(1.12) 
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where we used the identity 



1 — cos ( 



(1.13) 



to obtain Eq. ( 1.12| ) Transformation of 7i to real s pac e 
using the low q_L form of A'^ produces H^^ of Eq. ( ^.l[ ) . 
The absence of any qy dependence in the coefficient of 
in A''(q ) can be traced the equality B'^^/B^id — ^B^z 
[Eq. ( 2.20| )] that results from integrating out /i" from 
the original model. Effective Hamiltonians for /i" (r) and 
Uy{r) can be obtained by respectively integrating out Uy 
and u" from 7i'°* and u" from H.^^ . The results are 



2 J (27r)3 



where 



A;:(q)=A,(q) 



|A.(q)P |A,(q) 



A''(q)=A,(q)- 



A,(q) AyicO 
|A,(q)P 

A^(q)-|A.(q)PM.(q)' 



(1.14) 
(1.15) 

(1.16) 
(1.17) 



In the hmit qj^ ^ with (B2da^ /Kv) q^, both of 

these Hamiltonian reduce to Hh in Eq. ( 3.3S 



APPENDIX B: CALCULATION OF («)^) 

In this Appendix, the expression for the sliding colum- 
nar displacement fluctuations given in (3.19) is derived. 
To do this, we evaluate the integral of the the Fourier 
transformed SC correlator Gzz(q) over all (/-space: 



(2^)3 Bq^,+Kqi + KyqXp{qya)- 



(2.1) 



where p{u) = 2(1 — cosu)/v?. The fluctuations diverge 
at small wavenumbers q ~ '^/L, where L is the system 
size. To calculate how the fluctuations scale with La,, we 



set L, 



oo and L,, 



Note that the SC form for 



Gzz{^) is valid only when 3> x* , where x* — a/^y 
and /iy = Ky/K. The flrst step in the calculation is 
to perform the integration over with — s- oo and 
then use the fact that Gzz{^) is an even function of q 
so that the remaining integrals run over only positive q^ 
and qy. Note that taking the A^ — » cxd limit docs not 
alter qx^y — > divergences. The resulting expression 



T 



dqx 

9= 



— / dqy 



1 



fJ-lqyPiqya) 



where Aj, = 7r/a, is made dimcnsionless by changing vari- 
ables to V = qxX* and w ~ qya. The In^ divergence 



of the displacement fluctuations can be seen immediately 
by looking at the qx,y —> limit of this expression. 

The SC displacement fluctuations can be written as 
the sum of a continuum term Ic that does not depend on 
p{w) and discrete term Id that does depend on p{w). 



{{u 



n\2\ 



2t:^ y/BK^ 
where Ic and Id are dcflned by 



[Ic + Id] 



(2.2) 



Ir. = 



dv 



dw 



dw 



We first focus on the continuum contribution 7^. The 
integral over w is straightforward; 



Ic= / 

EE /(I) - 



0??; 



[f{Aya,v)-f{Ly'a,v)] 



(2.3) 



7-(2) 



where 



f{x,v) = In 



' + \J x^ ^ 



(2.4) 



The integral over v in l'^'' can be evaluated by separating 
the function 



JiKjCL, v) = ln[2Aya] -I- In 



^ + ^V^l + (VA,a)^ 



2 2 

into a constant term and a term that is well-behaved at 



small V. We then insert this expression into (2.3) and 
find that 



/(I) = ln[2Aya] ln[A,L,] 



(2.5) 



x'L-^ V 



■In 



The flrst term diverges with system size Lx , and the sec- 
ond term is nondivergent. f{aLy^,v) can also be sepa- 
rated into a constant term and a term that depends on 
V. We then integrate f{aL~^,v) over v to flnd 



iP = ln[2aL-i] ln[A,L,] 



(2.6) 



X* L 



dv 

— In 

-1 V 




The Ly dependence in the integrand of the second term 
can be moved to limits of the integral by changing vari- 
ables to s = Lyv/a. In contrast to the previous expres- 
sion for /i^-* in (2.5), the large s part of the integral in 



q) diverges with system size. The divergence can be 
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isolated by adding and subtracting ln[s/2]/s. The result- 
ing expression, 



1 



) = ln[2aL-i] ln[A,Lj + - In 



— ( In 



1 1 /— 

2 + 2^ 



--In 

ln[s/2] 



has two terms that diverge and two terms that do not 
diverge with system size. Note that Ly/L^Hy is 0(1) 
since ~ Ly. We then subtract ji^-* from Ic^\ drop the 
nondivergent terms, and set Ly ~ to obtain 



Ic = ln[A^L^] ln[A;,i,y] - ^ In^ 



AxLy 



2^y 



= i ln^[2/ij^Aj,ia; 



for the continuum contribution to the displacement fluc- 
tuations. 

The discrete contribution is obtained by evaluating 



Fiv), 



where 



F{v) = 



dw 



(2.7) 



(2.8) 



In the definition of F(v), we can take the lower limit to 
zero since the small w part of integral is well-behaved. 
In contrast to the continuum term, the discrete term di- 
verges logarithmically with system size. To see this, we 
expand F{v) around v — and find F{v) « F{0) + av'^ + 
0{v'^), where F{0) — ln[4/7r] and a is a constant. Thus, 
Id = ln[4/7r] ln[A2,i^] plus terms that do not diverge with 
system size. To make the argument of the In term in Id 
match the In^ term in Ic, we add and subtract the con- 
stant ln[4/7r] ln[2^yAy/Ax] to find 

/d = ln[4/7r]ln[2/ij,Aj,L^]. (2.9) 

We then add Ic and Id and find the following expression 
for the displacement fiuctuations in the limit oo 
and ^ Ly-. 



{{u 



n\2\ 



ll In^ [^Lx/x*] 



(2.10) 



where Sn{t) was defined previously in ( [3.23 ) and 

dt r , 1 - cosnw 



-^niAxi Aj 



du- 



t Jo y/W+U^p(u) 



(3.2) 



X 1 arctan 



A^z* 



Note that the integral defining An{Ax, Az) does not have 
an infrared divergence as t ^ [L^ oo), and thus it 
is a well-defined number. In the final step, we isolate the 
In divergence in the first term of (3.1) to find 



g^^\na)^2llSniO)ln[An{Ax,A,)Lx/x*], (3.3) 
where 

A„(A„A,) =exp[cii) -fci^) (A,) -fc(f) (A,, A,)] (3.4) 

depends on the layer index n and the ultraviolet cutoffs 
with 



Snit) 



ci'^A^A.) 



A, a;* 



^„(0) 

dt Snjt) 
1 t Sn{0) 

An{Ax,A,) 
Sn{0) ■ 



(3.5) 
(3.6) 
(3.7) 



Note that An{Ax, A^) is well-defined in the limit A^; 
oo. 



APPENDIX D: CALCULATION g'^^r) 

In this Appendix, we evaluate the SC position correla- 
tion function 



(4.1) 



between two DNA strands located in layer n = and 
separated by r in the xz plane. For general separations, 
.g^^-'(r) cannot be expressed in closed form. The aim of 
this Appendix is to calculate ^''^''(r) along the special 
directions z = 0, x x* and x = 0, z ^ z*. 



APPENDIX C: CALCULATION OF g^'^>{na) 



1. Large x, Small z Limit 



In this appendix, we calculate the diverge nt p art of 
the position correlation function g'^^^(na) in ( [3.22| ). We 
perform the integration, switch to dimensionless vari- 
ables, and find 



g^'Hna)^2ll 



dt 



X* L 



^ jSnit) + An{Ax,A,) 



(3.1) 



The following exp ressi on for g'-^'>{x,0) is obtained by 
setting z to zero in ( |3.16 ): 



?(2)(x,0)=T 



d^q 



I ~ cos{qxx) 



(2^)3 Sg2 + Kq^x + Kyqlqlpiqya) 



(4.2) 
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The first step in the derivation of g^'^^x, 0) is to perform 
the integration over with ^ oo. The Qz integration 
yields 



(4.3) 



where 



T( A ^ f'^'^^ 1-C0s(q,x) r/^' 
I{x,Ax) = / dqx- 



dqy 



<& + QyPiQyX*) 



We then decompose I{x, Ax) = Ic{x, Ax) + Id{x, Ax) into 
continuum and discrete contributions as we did previ- 
ously in Appendix M, where 



Ic{x,Ax) 

Id{x,Ax) 



, l-COs{qxx) r/"* dqy 
dqx- ' 







l-cos{qxx) r/^* 
dqx / dqy 



qx 







Since the A^^ — s- oo limit is well-defined, we calculate 
/c(x) = /c(x,oo) and /^(x) = Id{x, co) and drop terms 
that depend on the finite ultraviolet cutoff. 

To calculate the continuum contribution, we first set 
qx ~ uqy and then v ~ qyX. These changes of variables 
yield 



where 



Ic{x) 



K{v) 



-K{v), 



°° 1 — cosfww) 
du- 



10 uvT"+ 

The strategy for calculating the In^ x term in Ic{x) is to 
isolate the part of K{v) that scales as In-y for large v. To 
this end, we write 



(4.4) 



(4.5) 



K{v) 



— [1 — cos(mw)] 



/ — [1 — cos('uw)] 
Jo u 



dul — cos{uv) 



(4.6) 



It is obvious that only the first term has the correct scal- 
ing; the remaining terms in K(v) are then separated into 
constants and functions of v that are well-behaved either 
as ?; — > or u — > cx). This partitioning leads to 



K{v) = \n{Bv) + K{v), 
where B — 2e'', 7 is Euler's constant, and 



(4.7) 



K{v) 



du 



cosu 



du 



u 

cos(uw) 



du cos[uv) 



u Vl + u^ 



(4.8) 



scales as l/v^ for large v. 

We then insert K{v) into (4^) and break the integral 
over V into small- and large-w parts to obtain 



dv r'^''^ dv 

I,{x)^ / —K{v)+ / —\n\Bl^ (4.9) 

Jo V Ji V 

7Tx/x* 1 

^"-X(.). 



Next, we evaluate the integral over v in the second term, 
collect constants, and find 



Icix) ^Un^[2e^Trx/x*]+Ax, 



(4.10) 



where 



Ax = -lln'[2e^]+ f-K{v)+ H -K{v). (4.11) 
2 Jo V Ji V 

The second and third terms in Ax are finite since K{v) 
scales as for small v in the former and there is phase 
cancellation from the cos(ww) factor at large v in the later. 

We will now calculate the discrete contribution to 
f/^^-* (2;, 0). The first step is to rewrite Id{x) in dimcn- 
sionless form: 



Id{x) 



dv 

— [l-cos{vx/x*)]F{v), (4.12) 



where F{v) was defined previously in ( p.8[ ). We next 
break the integral over v into small- and large-w parts 
and take the x ^ x* and Aj, — s- cx) limits to obtain 

Id{x)^F{Q) r —\l-cos{vx/x*)]+ r —F{v) 
Jo V Ji V 

''^ dv 



-[Fiv) ^ F{0)]. 



(4.13) 



Note that taking the x x* limit removed the 

cos{vx/x*) terms from the last two terms in (4.13) due 
to phase cancellations. It is again obvious that the first 
term scales logarithmically with x/x*, and thus 



Idix)^ ln[4/TT] In 



where 



Bx^ r-[Fiv)-Fm 
Jo " 



Bx, 



dv 



(4.14) 



F{v) (4.15) 



is a constant. The last step in the calculation of g{x,0) 
is to add the continuum and discrete terms, Ic{x) and 
Id{x). The final result is 
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Cx 



(4.16) 



where = 2{Ax + B^) - ln^[4/7r] - 2 ln[4/7r] ln[27r]. 
and B x are com puted later in this appendix, see Eqs. 
(4.34) and ( 1.35| ). By us ing t hese equations, we eventu- 
ally find that Cx in Eq. ( 4.16 ) vanishes; Cx = 0. 



2. Large z, Small x Limit 

The calculation of g*-^^ (0, z) is similar to the calcula- 
tion of 5*^^-' (a;, 0). The e xpres sion for z) is obtained 
by setting x to zero in ( 3.16 ). The first step in the calcu- 
lation is to perform the integration over with Kz ^ oo 
which yields 



g(2)(o,z) 



T 



271^ y^BK,, 



(4.17) 



where 



I{z,Ax) = 



dqx 
Ix Jo 



dqy- 



^ ql + qlp{qyX*) 



In what follows, we set A^, — > oo, drop terms that depend 
on the finite ultraviolet cutoff, and define I{z) = I{z, oo). 
The second step is to change variables to u = qx/qy and 
V = Xzqy and decompose I{z) = Ic{z) + Id{z) into con- 
tinuum and discrete terms, where 



1 r dv 
2 



du- 



1 



-VUy/l-\-u'-^ 



V JO 

dv 



Ic{z 

Id{z)^ I —[F{v,0)-Fiv,vz/z*)] 
Jo " 

with T = TT^z/z*, z* — a^/fiyX, 
du 



(4.18) 
(4.19) 



F{v,t) 



•\/f^ -I- v?-y(u) 



(4.20) 



and i^(w,0) is equivalent to Fiv) defined in Eq. (2^ 

We first focus on the continuum contribution to 
g^^''(2:,0). The integral over v in lc(z) can be broken 
into small- and large-w parts. 



Ic(z) - - / —J{V) + - / —J{V), 

2 Jo V 2 Ji V 



'^'^Z^' dv 



where 



Jiv) 



oo 

du — 



(4.21) 



(4.22) 



The strategy is to extract the part of J{v) that scales 
as In?; for large v. If J{v) ^ \nv for large v, Ic{z) will 



scale as ln^[z/z*] as expected. Note that J(v ) scales as 
for small v, and thus the first term in ( 4.21 ) is a finite 
constant. After some algebra, we find 



J{v) = ln[i:>w] + J{v) 
where D — 2e^ ^ 



(4.23) 



J[v) 



du — 

^/2v U Ji u 

1 



du- 



Uy/l + lP 



and 



$(m) 



1 -I- 2u^ 



(4.24) 



(4.25) 



We ca n now insert the expression for J{v) in (4.23) into 
(4.21) and obtain the continuum contribution 



Uz) = \ln' 



2e^7r2- 

z* 



A, 



(4.26) 



where 



^. = -iln^[2e1 + i f'^Jiv) 
4 2 Jo u 



1 

2. /I 



dv 



-Jiv) 



(4.27) 



is a constant. Note that J{v) decays exponentially for 
large v, and thus the third term in A^ is finite. 

We now concentrate on the discrete contribution to 
^^^^(z, 0). The integral over v in /^(z) = Z^' -I- ij^^ can 
also be broken into small- and large-w parts, where 



r(i) 



1 dv 



[F{v,0) ~ F{v,vz/z*)] 



(4.28) 



(2) . . . 

and is an identical expression except the limits on 
the integral over v run from one to infinity. To isolate 
the Inz term in we change variables to i = vz/z* 
and take the z '3> z* limit. In the large z limit, ( 4.28| ) 
becomes 



ij^^ ^ F(0, 0) In 



" Z ' 








+ Bz,^ 


-L 


-Z* . 





1 dv 



[F{v,0) ~ F{0,0)] 
(4.29) 



where F(0,0) = ln[4/7r] and 

Bz= r^[F(0,0)-nO,t)] 
is a constant. The large-u contribution to 

r ^ 

li V 



I. 



(2) 



-Fiv,0), 



jFiO,t) (4.30) 



(4.31) 
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is simply a constant when z ^ z* . We then collect the 
discrete contributions and find 



= In 





r(2) 




'A 


In 


" Z ' 










.Z* . 



(4.32) 



where = 3^+ with B^ given by Eq. ( |4.15| ) . 

The last step in the calculation of 5(^^(0,2) is to add 
Ic{z) and Id{z). The final result is 



.(2)(0,z) = /2(_in^ 



32e^ 



(4.33) 



where = 2(^2 + S^) - 2 ln2[4/7r] - 2 ln[4/7r] InpeTTr^]. 

an d a re ca lculated later in this appendix [See Eqs. 
(4.36) and (4.37)]. By using these equations, we find 



3. Calculation of A^, B^, A^, and 

In Sees. 1 and 2 of this Appendix, we anticipated that 
the numerical constants A,j:,Bx,Az, and B^ have the val- 
ues: 



A, 
Bx 

A. 



TT 

24' 

-^+41n2(2) 
-[ln(2)][ln(^)]- 

48 ' 



ln2(7r) 



— +61n2(2) + 27ln(2) 
-[ln(2)][ln(7r)] - 7ln(7r) - ln2(7r). 



(4.34) 

(4.35) 
(4.36) 

(4.37) 



In t his se ction we outline calculations yielding Eqs. ( 4.34 ) 
to (4.37). We begin by deriving more explicit formulas 
f or t h ese consta nts. Thus, for A^, we obtain, by Eqs. 



(^, (|4^), and (yi 

A-x - 

with 



\ ln2(2eT) + + + 743), (4.38) 



A^^) = - f dvHv)^—^ 
Jo 

poo 



V 

cos(u) 



dv ln(u) 
1 V 



(4.39) 



du 



" 
°° du 



ln(u) 



1 



ln(M) 



(4.40) 



and 



du 
u 



1 



du 1 



- 1 



(4.41) 



Likew is e, for the constant A^, we obtain, by Eqs. (4.22) 
d^H), (|425| ), and (|l27| ). 



i[4^) + iln^(2) + 2ln(2)] 
l[Af+jAi^)], 



(4.42) 



where 



dx ln{x 



1 - e 



dx\a.(x)- 



(4.43) 



42)= /'dwl^jMin(uyTT^) 



du 



u^ 
1 



■.\n{uy/l + u^), (4.44) 



UV 1 + 



u 



VT 



(4.45) 



with $(u) defined in Eg. (4.25 ). W e proceed to compute 
t he in tegrals in Eqs. ( |4.38| ) to ( 4.45[ ). Integrals ( ^.4l|) and 
(4.45) can be done by elementary integration methods 
yielding 



4^)=ln(2), 4^) = ^. 



(4.46) 



Other integrals are not elementary, and we outline their 
calculation in the following. Thus, we compute the inte- 
grals ^i^-* and ^i^-* first by showing that 



AW = ^(1) - — 
and, next, by showing that 
All) 



1 ffr 

2 [d^ 



(4.47) 



(4.48) 



where r(z) signifies the gamma function. 

To obtain Eq. (4.47), consider the complex function 
f{z) = e^^ ln(z)/2:. In the complex z = x+iy plane, f{z) 
is analytic inside the contour C made of the following 
four segments: Ci : [z — x;e < x < R], C2 ■ [z = 
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i?e'^;0 < < 7r/2], C3 : [z = iy; R > y > 0], and 
C4 : [z = ee'^;iT/2 > 9 > 0]. By applying the Cauchy 
residue theorem to the above /(z) along the contour C = 
Ci + C2 + C3 + C4 in the limit e ^ and R 00, we 



directly obtain the relation (|4.47l) between A^^^ and Az 
Next, we demonstrate Eq. ( 4.48 ) by differentiating the 
standard integral representation of the gamma function. 
This yields 



dx ln(a;)e 



dx ln^(a;)e 



(4.49) 



(4.50) 



Wc then rewrite Eq. (4.50) as 



I dz2 



+ 



dxie-"" ~ l)ln^(x) 



dx\T?{x)+ / dxe-''\T?{x) (4.51) 



and integrate by parts the first integral [by writing (e~^ — 
l)dx = d{l — — x),etc.] as well as the last integral [by 
writing e~^dx = d{—e~^), etc.]. Aft er few elementary 
integrations, this yields our Eq. (4.48) which, combined 
with Eq. (4.47), yields values of the integrals A^x^ and 
Ai^\ To complete this calculation we need also the value 
of the second derivative of r(z) a.t z — 1 that enters Eq. 
(4.48). To compute it, we use a relation from the theory 
of the gamma function: 



dr{z) 

dz 



r(z)*(z). 



(4.52) 



with 



vI/(z) = -7+(l--) 
^ z' 

^3 z + 2' 



z + 1- 



(4.53) 



r(l) = 1, {dT/dz)z=i = *(1) = -7. By differentiating 
Eq. (p2|). 



dz^ 



dV 

dz 



= 7 



z=l 

dz 



dz 



z=l 

(4.54) 



The use of Eq. ( 4.53| ) to compute d^{z)/dz at z — 1 
yields 



d^r\ _ 2 1 1 1 1 

d^L_,^^^ ^22 + 32 + 42 



= 7' + y (4.55) 



By Eqs. (|^), ( f4.47D , and ( [4.48D , we finally obtain 



9 9 9 9 



24' 



12 



(4.56) 



(2) 

Nex t, we proceed to compute the integral Ax in Eq. 
(4.40). For this purpose, consider the integrals 



1 



ln(u) 



du 



t u Vl + u- 



: ln(u) 



A^Ht)- f- 

Jo u 



- 1 



du 



For t = 1, A''x\t = 1) = A'^' is the desired integral 



t u Vl + u2 

(2) 



(4.57) 



(4.58) 



whereas A],'{t = 1) = A],' = ln(2), by ( [4.46[ ). Further, 
by Eqs. (|4.57|) and ([4.58|), dA^x\t)/dt = - \n{t)/t, and 



dA'i'\t) / dt = —1/t. By integrating these relations over 



(4.59) 
(4.60) 



for any t > 0. By Eqs. ( [4.571) to ( [4.60D , with t = e ^ 0, 



^(2) = lim 
= lim 



du 



1 



u Vl + V? 
°° rfu 1 



ln(u) + 
ln(e) 



u Vl + V? 

To proceed, it is useful to form the difference 



Using here Eqs. ( 4.61 ) and ( [4.62 ), we find 



A = / du — , / — 

/o uVl + u'^ Jo V 



1 - 



(4.61) 
(4.62) 

(4.63) 

(4.64) 



The integral over v here can be done by changi ng v t 
with t = ln[(Vl + v'^ + l)/2]. This reduces Eq. ([4.64[) to 



A = 



du- 



:ln[ 



VT+v^ + 1, 



(4.65) 



By changing variables via u t with t ~ hl[{^/l + + 
l)/2], we eventually obtain A in a form involving a well 
known integral. 



dt- 



1 12 



(4.66) 
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By Eqs. ( ^^6^ ), ( |4.63D , and ( [4.46D , we obtain 



12 



In' (2) 



(4.67) 



Equations (H.38|) , (|4.46|) , (H.56|) , and ( [4.67| ) yield our final 
result for the constant anticipated in Eq. (4.34). 

Next, we compute the integral Ai in Eq. (4.44). For 
this purpose, consider the integral 



Af\t) 



, 1 - $(u) , , r 

du — ln(MV 1 + u^) 

uvl ■ 



du- 



1 



: ln(w vl 



(4.68) 



t U\/l + V? 

For t = 1, Af\t = 1) = Af^ is the desired integral. It 
i s eas y to s how along the lines we used to derive Eqs. 



( [4.59D and ([4.60D , that 



for any i > 0. Thus, by Eqs. ( [4.68D and ( |4.69| ), with 

t = e^O, 



= hm 



1 



ln(uVl + u2) + i^ 



ln'(2) 



2 

(4.70) 



Next, by Eqs. ([4.70D and ([4.6lD , 



1 



:ln(v/l + u2)- i^^. (4.71) 



By the change of variables u t, with u = Ve^* — 1, we 
S*! + S'2 - ln2(2)/8, where Si and fi-s 



find A 



(2) 



(2) 



are two well known integrals. Si = dt t/{e* + 1) 
^2/12, and ^2 = /o°° dt t/{e^*-l) = \ \^ dxx/{e^-l) = 
i(7rV6). Thus, 



Af^ - ^ 
By Eqs. (|7|) and ( ^J7| ) 



TT 



ln'(2) 



^(2) = 

24 



57r2 31n2(2) 



(4.72) 



(4.73) 



Equations (|4.42|) , ([4.46|) , (|4.56|) , and ( [4.73| ) yield our final 
result for the constant A^ anticipated in Eq. ( 4.361 ). 

We now proceed to discuss calculations yielding the 
value s of t he nu meric al constants and quoted in 
Eqs. ( |I35| ) and ( |43^ ). For B^, by ( jilsl ) and (U), we 
obtain 



S.r = / dw 

JO 



(4.74) 



By using ( p^ ), i.e., A^i\t) = Ai^^-ln(t) = ln(2)-ln(0, 
we find 



Br 



dw 



ln(2) + ln(?i7^p(u;)) 



dw 



Wy/p{w) 

ln(2) +ln(w)" 
w 



(4.75) 



Next, in the first integral above, we change from w to z = 
wyjp{w) = 2sin(it;/2), whereas in the second integral 
we change the variable w into z. After rearranging the 
expression thus obtained, we find 



Bcr 



dz 



1 



" _^_ ln(2) + ln(z) 



ln(2) + ln(z) 



(4.76) 



For convenience, we change variables via z x — z/2 
and thus obtain 



B^ = B^^^ +21n(2)B 
-[ln(2)][ln(7r)] 



(2) 



with 



and 



ij(i) 





1 


/ dx 




h 


Wl-x^ \ 



InV) 
2 

31n2(2) 



ln(a;). 





1 


/ dx 




'0 





(4.77) 



(4.78) 



(4.79) 



A similar formula can be derived for the constant B, 



Bz + Bx i ntrod uced in Sec. |D 2| of this Appendix. By 
Eqs. ( p^ and (p3o|)), we find for B^ 



B, 



^0 

(4.80) 



By treating Eq. (1.8C) in exactly the same way we treated 
above Eq. ( 4.75| ), we eventually find 



B^ = + [ln(2) +7]B 



(2) 



In' (2) 



-7ln(2) 



ln2(7r) 



7 ln(7r). 



(4.81) 
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To complete our calcu lation, we n eed the int egrals B^^^ 
and in Eqs. (^J8|) and (^J9|) . B'-^\ Eq. can 
be calculated by elementary integration methods yielding 



5(2) = in(2). 



(4.82) 



On the other hand, the integral B^^\ Eq. (4.78), is not 
elementary. As detailed below, we find 



7r2 In" (2) 
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(4.83) 



The value s of and Bz = Bz + B^ quote d in Eq s. (|4.35| ) 
and (H.37|) , directly follow from Eqs. ( |4.77| ), ( [4.81D , ( p^ , 
and ( ^4.83[) . It remains to outline the calculation yielding 
the value of B'^' in Eq. (4.82). It is obtained by relating 
t he in tegral B^'^\ Eq. ( [4 781) to the integral A^^\ Eq. 



( [4.40D , computed in Eq. ( ^.67|) . We find that 



(4.84) 



To derive relation (4.84), we apply the Cauchy residue 
theorem to the complex function f[z) — ln(z)/z\/l — z'^ 
along the sam e contour that has been used before to de- 
rive Eq. ( [4.47| ) [see the text following Eq. (f4.48| )]. This 
yields the relation 



dx i_ 

Je X JT^ 



a/T-cc 

^2 /-oo 



: ln(x) 



8 



du 



:ln(u) + 0(e), 



(4.85) 



where 0(e) ^ as e ^ 0. Equation (4.85) is identical to 



dx 

X 



ln(x) 



du 



ln(u) 



+ 0(e). (4., 



By r ecalli ng here Eq. ( 4.61 ) and taking the lim it e in 
Eq. ( 4.86 ), we eventually obtain the relation ( [4.84 ) used 
to compute B^-^\ 



APPENDIX E: CALCULATION OF g^^^r,na) 

In this appendix, we will outline the evaluation of Eq. 
( 3.25 ) for g^'^^(r, na) = g^^^^ {x , z , no) . We begin with the 
expression, Eq. (3.17) for g'^^ (r, na). We assume the con- 
tinuum limit Az 



00 and integrate over to obtain 



5(3)(r,na) = 2/2 



dt 



du- 



1 — cos(r7,M) 



t Jo ty/t^ +u'^p{u) 
1 - cos(te/a;*)e-(*^/^*^V*'"'+P(")l . (5.1) 



This expression can now be used to evaluate various lim- 
its. 



1. a; ^ a-* and z — 
From Eq. (|5.lD, we have 

5(3)(x,0,na) = 2ll / ^3^1 - cos(te/a:*)], 
Jo r 



(5.2) 



where Sn{t), defined in Eq. ( 3.23 ), i s pr oportional to 1/t 
for i ^ 1. Thus the integral in Eq. (5.2) is convergent at 
large and we can take the continuum limit A^^ — > 00. 
Thus, we have 

r°° df 

g^'Hx,Q,na) = 2ll / 5„(t)[l - cos(te/x*)] 

2ll / '^[l-COs{tx/x*)]Sn{t) 







t 



dt, 



^211 J j[l-cos{tx/x*)]Snit) (5.3) 
2ll[ln{x/x*)+j]S,M 

df 

+ 2ll / -[l^COs{tx/x*)][Sn{t)~Sr.{Q)] 



'21. 



dt 



[l~COB{tx/x*)]Sn{t) (5.4) 



The contributions from the cos[tx / x*) terms in the inte- 
grals in Eq. ([3.4|) vanish when x/x* 3> 1, leaving 



5(3)(a;,o,na) = 2llS^{0)\n{D,,x/x*), 



where 



(5.5) 



(5.6) 



with Cn'^ and Cn\^x) given, respectively, by Eqs. ( |3.5| ) 
and (3.6). 



2. a; = 0, 2: > z* 

Again, we will simplify our discussion by ta king the 
limit Aa; — > 00 and integrating over in Eq. ( p.25| ) to 
obtain 



g(3)(0,z,na) = 2ll / dt I du 



00 clX 



1 — cos(nu) 



t^t"^ + u^-p{u) 



1 - e 



-{tzlz')^t^+v?y{u) 



(5.7) 



Changing variables via t = u^piu)., we obtain 



.(3)(0,z,na) = 2ll ^-^^J[{z/z*)uMu)], 



y/u'^p{u) 



(5.8) 
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where J{v) is defined in Eq. (4.22). Next, using J(v) 
\n{2e'') + J{v) [Eq. (|4.23|)], we obtain 



where 



g(3)(0,z,na) = 2?2[ln(z/z*) +7 + ln2] / du 



^ 1 — cos(nM) 



$(ri,r) = ig2^(u^(r,ri)+{i^(ri,0) + aa,ui(ri)2) 



= 5^2)(r-ri,0) + .g^^^(ri,0) 
+.g(2) (r, a) + 5(2) (Q, a) - 5(2) _ ^ , a) - 5(2) (r, , „) 

+a(a.«i(ri)[{ii(r,ri)-^i°(ri,0)]); (6.4) 



^ u2p(u) 



2 , 1 - cos(nu) , . 2 ^ M 
+2i„ / du , „ , . ln[u p[u)\ 



+211 rdu^^^°'("") ^-./.*^.,2 



J[(z/z*)u^p(u)] (5.9) 



In the hmit z/z* —s- cx), J[(z/z*)u^p(u)] tends to zero for 
ah u, and we obtain 



7(3)(0,z,na) = 2Z2s„(0) ln(i;„z/z*), 



(5.10) 



where 

Eu — 2e^ exp 



1 - cos(nu) , r 2 / M 



Sn{Q) Jo v/^^M^ 



(5.11) 



APPENDIX F: INTER-PLANE CORRELATIONS 
WHEN / 



In this appen dix, we will derive Eq. ( 3.4(: ) from the 
expression, Eq. (3.41) for 5(r,na). We begin by looking 
at the linear term in in 5(r, a): 



S{r,a) ^ V^"/2T J dVi(e~*'°"-(°) 

X ^gi'lo[uliri)-ul{ri)]^iqoaa^ul{ri) ^ ^^^^^ ^iqoulir)^^ ^g^-j 



Xg«9o[u°(ri)-«°(0)]gigoa9^«]^(ri)^^^ 



2T 



(f2^^^gigo[«^(r)+«i(ri)] 



(6.2) 



where the subscript SC indicates that the averages are to 
be evaluated with respect to the sliding columnar Hamil- 
tonian, H^'^ of Eq. ( pl9| ). Since Ti^'^ is quadratic in 
u" (r) and Uy (r) , the averages in this equation can be per- 
formed exactly, and 5*2 (r, a) can be expressed as an ex- 
ponential of correlations functions of M"(r) and dzUy{r). 
The second term in Eq. ( |6.2[ ) has terms in the exponen- 
tial proportional to — qo(['u"(r)]^) for n = 0, 1, which di- 
verge in thermodynamic limit and cause the exponential 
to vanish. Thus only the first term of Eq. (6.2) survives, 
and we have 



(6.3) 



here, g^^'(r, na) is defined in Eq. ( |3.1S| ), and 

u^(r,r') = <(r)-<(r'). (6.5) 

This expression is quite complex. It, however, simpli- 
fies considerably if we set Ky = and But = 0. Then 
^'^■'(r, a) = 0, and all cross terms in u"(r) and u™(r') 
vanish. In this case. 



<firur)^g^^Hr-r^)+g^^^ir^) + Wy 

,(2)( 



(6.6) 



where g^^^r) = g'^^\r,0) [See Eq. ( p^ j and Wy 



The generalization of Eqs. ( |6.3D and (3.4) to n > 1 is 
straightforward. The leading contribution to ^(r, no) is 



5(r, no) 
where 



yu 

2T 



J d2ri...d2r„e-*('-i' -''--'^\ (6.7) 



\ \m=0 m=l / I SC 

(6.8) 

where ro = and r„+i = r. This function can be 
expressed in terms of the reduced correlation functions 
(7'-^''(r,n) for u"(r) and correlation functions involving 
dzU^iy). In general, it will have terms connecting all 
pairs of layers, and the evaluation of S{r, na) is highly 
nontrivial. 

We can obtain a useful approximation by setting Ky = 
to eliminate all g'^' (r, na) for n ^ and by setting 
Bull = to eliminate couplings between u" and u™ . Even 
in this approximation, couplings between distant layers 
arise from (SzU^ (r)9zU™(r')). This function, however, 
dies off rapidly with n—m, and we will set it equal to zero 
when n ^ m. [In another publication [ p^ , we will use 
a variational procedure to calculate S'(r,na). The form 
of 5(r,na) evaluated with this procedure is very similar 
to that obtained using the above approximations.] With 
these approximations, 

n 

$(ri,...,r„,r) ^ g^'^) ~ Tm+i) + nWy (6.9) 

m=0 

and 
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S{r,na) = ^ j d\^...d\,,S2{vi) 

x^2(ri-r2)...52(r„-r), (6.10) 

where 

S2{r) = (e**[<W-<Wl) (6.11) 
is the inplane density correlation function and where 

^" = ^"6-^". (6.12) 



Fourier transforming this equation, we obtain Eq. ( 3.4(]| ) 
in the text. 



c: 



dy 



dy^. 

1 y 



(7.7) 



and 



J, 



Fn{y) — -f- I cos(nu) -yit^^^u) exp ~y\/ u^p{u) 



n Jo 



Note that C^'^{Ax) diverge with A^^, and thus £'„(r) does 
not have a well-defined A^. ^ oo limit. 



APPENDIX G: INTERACTION ENERGY 
BETWEEN EDGE DISLOCATIONS 

In this appendix, we evaluate the interaction energy 
E{r, no) between two dislocations with separation x ~ 
(r,na) in the limits x ^ x* , z — Q and z ^ z* , x = Q. 
To obtain the x^ x* , z — Q limit, we must evaluate 



En{x) 



T 



X cos[riM] (1 — cos[tx/a;*]) . 



(7.1) 



We isolate the In x divergence by adding and subtracting 
v?p{u) under the u-integral. This procedure yields 



for x^ X* . In the above expression, {Ax 
withS„(A,) = A,x*e^"(A^), 



Bn{Ax) 



dt 



Jn{t) 



and 



Jn{t) = / duy^t'^ + u'^p{u) COS nu. 



(7.2) 

(7.3) 
(7.4) 



The z ^ z* , X = limit is obtained in a similar way 
from 



Eniz) 



dt 



t JO 

X cos[nu] (^1 ~ exp 



du\/ + vP'piu) 



:ty/t^ + u^p{u) . (7.5) 



We find that 

En{z) = -Jn In [Q(A,)|z|/z*] (7.6) 

scales logarithmically for z ^ z* , where C^{Ax) — 
S„(A,)e^"(A^), 
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